Periodic-Orbit Theory of Universality in Quantum Chaos 
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We argue semiclassically, on the basis of Gutzwiller's periodic-orbit theory, that full classical 
chaos is paralleled by quantum energy spectra with universal spectral statistics, in agreement with 
random-matrix theory. For dynamics from all three Wigner-Dyson symmetry classes, we calculate 
the small-time spectral form factor K (r) as power series in the time r. Each term t" of that series is 
provided by specific families of pairs of periodic orbits. The contributing pairs are classified in terms 
of close self-encounters in phase space. The frequency of occurrence of self-encounters is calculated 
by invoking ergodicity. Combinatorial rules for building pairs involve non-trivial properties of per- 
mutations. We show our series to be equivalent to perturbative implementations of the non-linear 
sigma models for the Wigner-Dyson ensembles of random matrices and for disordered systems; our 
families of orbit pairs are one-to-one with Feynman diagrams known from the sigma model. 

PACS numbers: 05.45.Mt, 03.65.Sq 



I. INTRODUCTION 
A. Background 

In the semiclassical limit, fully chaotic quantum sys- 
tems display universal properties. Universal behavior has 
been observed for many quantities of interest in such dif- 
ferent areas as mesoscopic transport or nuclear physics. 
One paradigmatic example stands out and will be the 
object of our investigation: According to the Bohigas- 
Giannoni-Schmit (BGS) conjecture put forward about 
two decades ago pj , highly excited energy levels of generic 
fully chaotic systems have universal spectral statistics. 
This conjecture is supported by broad experimental and 
numerical evidence |2|. L|. 

Level statistics can be characterized by the so-called 
spectral form factor. The level density p(E) = J2 i 5{E — 
Ei) of a bounded quantum system [Ei denoting the en- 
ergy levels) is split into a local average p(E) and an os- 
cillatory part posc(E) describing fluctuations around that 
average. The form factor is defined as the Fourier trans- 
form of the two-point correlator (p osc (E + | )p OS c(E— §)) 
w.r.t. the energy difference e, 



K(r) 



de 



E- 



(1) 



here the time r, conjugate to the energy difference, is 
measured in units of the so-called Heisenberg time 



T H = 2irhp(E) 



n(E) 



(2) 



with fl(E) denoting the volume of the energy shell and 
/ the number of freedoms. Since the study of high- 
lying states justifies the semiclassical limit, we may take 
h — ► Q,Th — > oo, for fixed r. To make K(t) a plottable 
function, two averages, (...) in lfl]l. are necessary, like 
over windows of the center energy E and a small time 
interval At <§; 1. 



Given full chaos, K(r) is found to have a universal 
form, as obtained by averaging over certain ensembles 
of random matrices [2, 0, H ■ In the absence of geomet- 
ric symmetries, the prediction of random-matrix theory 
(RMT) only depends on whether the system in question 
has no time-reversal (T) invariance (unitary case), or 
is T invariant with either T 2 = 1 (orthogonal case) or 
T 2 = -1 (symplectic case). RMT yields for < r < 1 



K{t) = < 



T 

2t 



- rm(l 
fln(l-r 



unitary; 
2r) = 2r - 2r 2 + 2r 3 - . . 

orthogonal; 



symplectic. 



(3) 

However, a proof of the faithfulness of individual 
chaotic dynamics to random-matrix theory, and even the 
assumptions required for a proof, have thus far remained 
a challenge. In the present paper, we take up the chal- 
lenge and derive the small-T expansion of K(r) for in- 
dividual systems; we employ ergodicity and hyperbolic- 
ity of the classical dynamics. Moreover, we require all 
classical relaxation times (related to Ruelle-Pollicott res- 
onances and Lyapunov exponents) to be finite; we need 
this property to make sure that even the shortest quan- 
tum time scale of relevance, the so-called Ehrenfest time 



Te ~ In c ™ is much larger than any classical time scale. 

Following H, 0, 0i we start from Gutzwiller's trace 
formula which expresses the level density as a sum 
over classical periodic orbits 7, 



IE) - -4 ^ V vV 1 ^, 



(4) 



wherein A 1 is the stability amplitude (including the 
Maslov phase) and S 1 the action of the 7th orbit. By 



Q, the form factor becomes a double sum over orbits, 
K(t) = ±/j2A,A;,ei^-^S^T H - (5) 

\7,7' / 

Tj is the period of 7. For h — > 0, only families of or- 
bit pairs with small action difference \S-y — Sy| ~ can 
give a systematic contribution to the form factor. For 
all others, the phase in (JSJ oscillates rapidly, and the 
contribution is killed by the averages indicated. Fluctu- 
ations in quantum spectra are thus related to classical 
correlations among orbit actions @- The first periodic- 
orbit approach to K(t) was taken by Berry p|, who de- 
rived the leading term in J3J using "diagonal" pairs of 
coinciding (7' = 7) and, for time T-invariant dynamics, 
mutually time-reversed (7' = T7) orbits, which obvi- 
ously are identical in action. Starting with Argaman et 
al. |(|, off-diagonal orbit pairs were studied in |g. lTollll| . 
The potential importance of close self-encounters in orbit 
pair s was first spelled out in work on electronic transport 
Il2| and qualitatively discussed for spectral fluctuations 

The family of orbit pairs responsible for the next-to- 
leading order was definitely identified in Sieber's and 
Richter's seminal papers Q for a homogeneously hyper- 
bolic system, the Hadamard-Gutzwiller model (geodesic 
motion on a tcssclatcd surface of negative curvature with 
genus 2). Their original formulation was based on small- 
angle self-crossings of periodic orbits in configuration 
space. In each pair, the partner 7' differs from 7 only 
by narrowly avoiding one of its many self-crossings. The 
two orbits almost coincide in one of the two parts sepa- 
rated by the crossing, while they are nearly time-reversed 
in the other part. In phase space, both orbits contain an 
"encounter" of two almost time-reversed orbit stretches. 
They differ only by their connections inside that en- 
counter; see Fig^ 

As shown in |14| . Sieber's and Richter's reasoning can 
be extended to general fully chaotic two- freedom systems. 
One partner orbit 7' arises for each encounter. The action 
difference within each orbit pair 0, 0, 0] can be de- 
rived using the geometry of the invariant manifolds |l7j . 
It thus turned out helpful to reformulate the treatment 
in terms of phase-space coordinates 0] , which may 
also be applied to systems with more than two freedoms 

E- 

In [jjj, we showed that the r 3 -contribution to the form 
factor originates from pairs of orbits which differ either 
in two encounters of the above kind, or in one encounter 
that involves three orbit stretches. 

In the present paper we demonstrate how the whole 
series expansion of K (r) is obtained from periodic orbits. 
Beyond furnishing details left out in our previous letter 
|20| we here cover systems with more than two freedoms 
and from all three Wigner-Dyson symmetry classes. For 
the symplectic case we employ ideas presented in 0, . 
For related work on quantum graphs, see for the first 
three orders of K(t) and [24| for a complete treatment. 




FIG. 1: Sketch of a Sieber/Richter pair in configuration space: 
The partner orbits, depicted by solid and dashed lines, dif- 
fer noticeably only inside an encounter of two orbit stretches 
(marked by antiparallel arrows, indicating the direction of mo- 
tion). The sketch greatly exaggerates the difference between 
the two partner orbits in the loops outside the encounter and 
depicts the loops disproportionally short; similar remarks ap- 
ply to all subsequent sketches of orbit pairs. 



B. Overview 

We set out to identify the families of orbit pairs re- 
sponsible for all orders of the r-expansion. The key 
point is that long orbits have a huge number of close 
self-encounters which may involve arbitrarily many or- 
bit stretches. We speak of an Z-encounter whenever I 
stretches of an orbit get "close" in phase space. "Close- 
ness" will be quantified below such that we may speak of 
the beginning, the end, and the duration of an encounter. 
Fig. |3i highlights two such encounters inside a periodic 
orbit, one 2-encounter and one 4-encounter. Here, as al- 
ways, we sketch orbit pairs in configuration space, with 
arrows =E indicating the direction of motion inside the 
encounter stretches. The relevant encounters will turn 
out to have durations of the order of the Ehrenfest time 
Te ~ In c ™ st - ; even though logarithmically divergent in 
the semiclassical limit (and thus larger than all classical 
time scales), these encounter durations are vanishingly 
small compared to the orbit periods, which are of the 
order of the Heisenberg time, T 7 ~ Th ~ h~f +1 . In 
between different self-encounters an orbit goes through 
"loops" , represented by thin full lines. 

Self-encounters are of interest since they lead us from a 
periodic orbit 7 to partners 7' which differ from 7 notice- 
ably only inside a set of encounters (see the dashed orbit 
in 01). In contrast, the orbit loops in between encoun- 
ters are almost identical. The almost coinciding loops of 
7 and 7' are differently connected inside the encounters. 

Not all reshufflings of connections inside an encounter 
yield a partner orbit. For example, reconnections as in 
Fig. |2t> give rise to a "pseudo-orbit" decomposing into 
three separate periodic orbits. Pseudo-orbits are not ad- 
mitted in the Gutzwiller trace formula. 

For T-invariant dynamics, we also must account for 
encounters whose stretches only get close up to time re- 
versal as in =± and =^; see Fig. 01 Correspondingly, 
loops inside mutual partner orbits may be related by time 
reversal. 

We thus obtain a natural extension of Berry's diagonal 
approximation. Instead of considering only pairs of orbits 
which exactly coincide (or are mutually time-reversed), 
we employ all pairs whose members are composed of sim- 
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FIG. 2: a) Solid line: Periodic orbit 7 with one 4-encounter 
and one 2-encounter highlighted by bold arrows. Dotted line: 
Partner 7' differing from 7 by connections in the encounters, 
b) Other reconnections yield a "pseudo-orbit" decomposing 
into three periodic orbits (dashed, dotted, and dash-dotted). 




FIG. 3: Periodic orbit 7 with one 2- and one 3-encounter 
highlighted, and a partner 7' obtained by reconnection; the 
encounters depicted only exist in T-invariant systems. 



ilar (up to time reversal) loops. 

We proceed to classify these orbit pairs. Partner or- 
bits may differ in a number vi of I -encounters; we shall 
assemble these numbers to a "vector" v — (1)2,1)3, .. .). 
The total number of encounters is given by V = X^>2 Vl - 
The number of orbit stretches involved in encounters, 
coinciding with the number of intervening loops, reads 

L = J2i>2 lv l- 

The orbit pairs (7,7') related to a fixed vector v may 
have various structures. Each structure corresponds to 
a different ordering (and, given T-invariance, different 
sense of traversal) of the loops of 7 inside the partner 
orbit 7'. When drawing orbit pairs as in Figs. [21 and [21 
these structures differ by the order in which encounters 
are visited in the original orbit 7, and by the relative di- 



rections of the stretches within each encounter (i.e. =E 
vs. 5=^ , or 3 vs. =± in time-reversal invariant sys- 
tems). Moreover, different reconnections inside the same 
encounter may give rise to different partners, and hence 
different structures. We will see that structures have a 
one-to-one correspondence to permutations, which will 
be used in Section IlIII to determine the number N(v) of 
structures related to the same v. 

Orbit pairs sharing the same v and the same structure 
may still differ in the phase-space separations between 
the encounter stretches. We shall parametrize those sep- 
arations by suitable variables s, u and determine their 
density wt{s,u) inside orbits of period T. The double 
sum JSJ over orbits defining the spectral form factor will 
be written as a sum of contributions from families of orbit 
pairs, with the family weight oc N (v)wt (s , u) . 

This paper is organized as follows. To free the presen- 
tation of unnecessary details we mostly disregard compli- 
cations due to / > 2 and "non-homogeneous" hyperbol- 
icity (i.e. Lyapunov exponents different for different peri- 
odic orbits). In Section [HI we will study the phase-space 
geometry of encounters and derive the density wt{s,u). 
The purely combinatorial task of determining the num- 
ber of structures N(v) is attacked in Section IIIII with 
the help of the theory of permutations. We thus obtain 
series expansions of K(t) for individual chaotic systems 
with and without T-invariance; those series fully coin- 
cide with the RMT predictions for the Gaussian orthog- 
onal ensemble (GOE) and the Gaussian unitary ensem- 
ble (GUE), respectively. Section ITVl generalizes these re- 
sults to systems where spin dynamics accompanies trans- 
lational motion; in particular, we find agreement with the 
Gaussian Symplcctic Ensemble (GSE) given T-invariance 
with T 2 = -1. In Section |V| we show that our semiclas- 
sical procedure bears a close analogy to quantum field 
theory. In fact, our families of orbit pairs are equivalent 
to Feynman diagrams met within the theory of disor- 
dered systems and the perturbative implementation of 
the so-called nonlinear sigma model. Finally, we present 
conclusions in Section IVII Further details, including a 
generalization to / > 2 and non-homogeneous hyperbol- 
icity, and remarks on the action correlation function of 
@, are given in Appendices. 



II. PHASE-SPACE GEOMETRY OF 
ENCOUNTERS 

A. Fully chaotic dynamics 

At issue are fully chaotic, i.e., hyperbolic and ergodic 
Hamiltonian flows without geometric symmetries with 
/ = 2 "classical" freedoms. In the orthogonal case, the 
Hamiltonian is assumed to be T invariant, THT^ 1 = H , 
with an anti-unitary time-reversal operator T squaring 
to unity. For convenience, we assume T to be the con- 
ventional time-reversal operator T(q, p) = (q, — p); that 
assumption does not restrict generality, since all Hamilto- 
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nians with non-conventional time-reversal invariance can 
be brought to conventionally time-reversal invariant form 
by a suitable canonical transformation [l7j . 

For each phase-space point x = (q, p) it is possible 
to define a Poincare surface of section V orthogonal to 
the trajectory passing through x. Assuming a Cartesian 
configuration space (and thus a Cartesian momentum 
space), V consists of all points x + Sx = (q + £q, p + 5p) 
in the same energy shell as x whose configuration-space 
displacement 5q is orthogonal to p. For / = 2, V is 
a 2-dimensional surface within the 3-dimensional energy 
shell. Given hyperbolicity, V is spanned by one stable di- 
rection e s (x) and one unstable direction e"(x) [2S|. We 
may thus decompose 5x as 

<5x = se' s (x)+we n (x). (6) 

As long as two trajectories passing respectively 
through x and x + 5x remain sufficiently close, we may 
follow their separation by linearizing the equations of mo- 
tion around one trajectory, 

s(t) = A(x,<)" 1 s(0) 

u(t) = A(x,<)m(0). (7) 

Here, s(t) and u(t) denote stable and unstable compo- 
nents in a co-moving Poincare section at x(i) , the image 
of x = x(0) under time evolution over time t. In the long- 
time limit, the fate of the stretching factor A(x, t) and 
thus of the stable and unstable components is governed 
by the (local) Lyapunov exponent A(x) > 

A(x,i)~e A W t . (8) 

The x-dependence of A and A will be relevant only in 
Appendix IB II when we treat non-homogeneous hyper- 
bolicity; until then, we may think of these quantities as 
constants. As in 0, El [Si the directions e s (x) and 
e" (x) are mutually normalized by fixing their symplectic 
product as 

e*(x) A e s (x) = e"(xf ( J ^ ) e * (x) = 1 ' (9) 

In ergodic systems, almost all trajectories fill the cor- 
responding energy shell uniformly. The time average of 
any observable along such a trajectory coincides with an 
energy-shell average. 

Periodic orbits are exceptional in the sense that they 
cannot visit the whole energy shell. However, long pe- 
riodic orbits still behave ergodically: According to the 
cquidistribution theorem [26| (see also Appendix IB II) . a 
time average over an orbit 7 augmented by an average 
over all 7 from a small time window, with the squared 
stability coefficient as a weight, equals the energy-shell 
average with the Liouville measure. A special case is the 
sum rule of Hannay and Ozorio de Almeida |27j 




Ergodicity makes, in the limit of long times, for a uni- 
form return probability: A trajectory starting at x again 
pierces through V in a time interval (t, t + dt) with stable 
and unstable components of x(t) — x(0) (or Tx(t) — x(0)) 
lying in intervals (s, s+ds) 7 (u, u+dii) with uniform prob- 
ability ^dsdudt. 

B. Encounters 

To parametrize an ^-encounter, we introduce a 
Poincare surface of section V transversal to the orbit at 
an arbitrary phase-space point Xi (passed at time t\) in- 
side one of the encounter stretches. The exact location of 
V inside the encounter is not important. The remaining 
stretches pierce through V at times tj (j — 2, . . . , I) in 
points x.j. If the j-th encounter stretch is close to the 
first one in phase space, we must have Xj w Xi; if it is al- 
most time-reversed with respect to the first one, we have 
Txj ps xi. In the sequel, we shall shorthand as yj pa xi 
with yj either Xj or Tx-,- . 

The small difference yj — xi can be decomposed in 
terms of the stable and unstable directions at xi , 

yj - Xi = Sj-e s (xi) + Uj-e u (xi) ; (11) 

the stable and unstable components Sj, Uj depend on 
the location Xi of the Poincare section V chosen within 
the encounter. If we shift V through the encounter, the 
stable components will asymptotically decrease and the 
unstable components will asymptotically increase with 
growing t\, according to 105}. 

We can now give a more precise definition of an l- 
encounter. To guarantee that all I stretches are mutually 
close, we demand the stable and unstable differences \§j\, 
\v,j\ of all stretches from the first one to be smaller than 
a constant c. The bound c must be chosen small enough 
for the motion around the I orbit stretches to allow for 
the mutually linearized treatment (JJJ; however, the exact 
value of c is irrelevant. 

The stable and unstable coordinates determine the du- 
ration t cnc of an encounter. We have to sum the dura- 
tions of the "head" of the encounter (i.e. the time t u until 
the end of the encounter, when the first of the unstable 
components \iij \ reaches c) and its "tail" (i.e. the time t s 
passed since the beginning of the encounter, when the last 
of the stable components \Sj\ has fallen below c). Using 
the exponential divergence of the unstable phase-space 
separations (J8J, we see that the coordinates \v,j\ approx- 
imately need the time j In t#t to reach c; similarly, the 

stable coordinates need 4- In t#t . We thus obtain 

. f 1 , c 1 . ri, c i 

t u = mm <^ - In — — ■ } , t s = mm <^ - In — — > , 
j I A \uj\) j [X \sj\) 

1 c 2 
tcnc t s + t u — In pn a TT j (12) 

A maxj{|Sj|} maXj{|Wj|} 

in view of Q , the duration t cnc remains invariant if the 
Poincare section V is shifted through the encounter. 
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FIG. 4: Connections between left and right ports in partner 
orbit 7'. In b), the encounter splits into three pieces respec- 
tively containing the two upper, middle, and lower stretches. 



An Z-encounter involves / different orbit stretches 
whose initial and final phase-space points will be re- 
ferred to as "entrance" and "exit ports". If all encounter 
stretches are (almost) parallel, as in all entrance 
ports are located on the same side of the encounter, and 
the exit ports are located on the opposite side. If the en- 
counter involves mutually time-reversed orbit stretches 
as E±, this is no longer the case. Thus, it is useful to 
introduce the following convention: All ports on the side 
where the first stretch begins are called "left ports" , while 
those on the opposite side are " right ports" . For paral- 
lel encounters, "entrance" and "left" are synonymous, as 
well as "exit" and "right" . 



C. Partner orbits 

The partner orbits (7,7') differ from one another only 
inside the encounters, by their connections between left 
and right ports. We shall number these ports in order of 
traversal by 7, such that the i-th encounter stretch of 7 
connects left port i to right port i. Inside 7', the left port 
i is connected to a different right port j; see Fig. 

We must reshuffle connections between all stretches of 
a given encounter. In contrast, Fig. 0Jd shows reconnec- 
tions only between stretches 1 and 2, stretches 3 and 4, 
and between stretches 5 and 6 of a 6-encounter which 
therefore decomposes into three 2-encounters. 




-^left port 2 



right port 2 




unstable 



FIG. 5: Piercings Xi, y2 of 7 (full line) and piercings y[, y' 2 
of 7' (dashed line) for a Sieber/Richter pair, depicted a) in 
configuration space, with arrows indicating the momentum 
of the above phase-space points, and b) in the Poincare sec- 
tion V parametrized by stable and unstable coordinates. The 
symplectic area of the rectangle is the action difference AS. 

stable 



S 3 

s 2 - 



*i iy/ 



unstable 



U 



1. Piercing points 

The partner 7' also pierces through our Poincare sec- 
tion V . The corresponding piercing points are deter- 
mined by those of 7. In particular, the unstable coor- 
dinates of a piercing point depend on the following right 
port. If two stretches of 7 and 7' lead to the same right 
port, they have to approach each other for a long time - at 
least until the end of the encounter (which has a duration 
~ Te) and half-way through the subsequent loop. Hence, 
their difference must be close to the stable manifold, and 
the unstable coordinates almost coincide. Similarly, the 
stable coordinates are determined by the previous left 
port, since stretches with the same left port approach for 
large negative times. If a stretch of 7' connects left port 
i to right port j, it thus pierces through our Poincare 



FIG. 6: Piercing points of 7,7' differing in a 3-encounter; 
inside 7', left ports 1,2,3 are connected to right ports 2,3,1, 
respectively. 



section with stable and unstable coordinates 

s'i ~ u'i ~ Uj. (13) 

For instance, if 7 and 7' differ in a 2-encounter, 7' 
connects left port 1 to right port 2, and left port 2 to 
right port 1; see Fig. |SJi. Thus, the encounter stretches 
of 7' pierce through V in with s[ ~ §1 — 0, u[ rs U2, 
and in y' 2 with s' 2 ~ §2, u' 2 ~ u\ = 0, which together with 
the piercings of 7 span a parallelogram in phase space 0J 
(a rectangle in Fig. [SJd, by artist's license). In Fig. HJ1 we 
visualize the locations of yj inside V for a 3-encounter. 
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2. Action difference 

We can now determine the difference between the ac- 
tions of the two partner orbits, first for 7, 7' only dif- 
fering in one 2-encounter. Generalizing the results for 
configuration-space crossings in we will show that 

the action difference is just the symplectic area of the 
rectangle in Fig. |Hp 0, E3|. Consider two segments 
of the encounter stretches in Fig. [SJi, leading from the 
first left port to the piercing point Xi of 7, and to the 
piercing point y[ of 7', respectively. Since the action 
variation brought about by a shift dq of the final coor- 
dinate is p • dq, the action difference between the two 
segments will be given by AS^ 1 ) = J*, 1 p ■ dq. The in- 

y 1 

tegration line may be chosen to lie in the Poincare sec- 
tion; then it coincides with the unstable axis. Repeat- 
ing the same reasoning for the remaining segments, we 
obtain the overall action difference AS = S 7 — Sy as 
the line integral AS* = <f p • dq along the contour of the 
parallelogram y[ — + Xi — > y' 2 — ► y 2 — * y^, spanned by 
y[ — xi = M2e tt (xi) and y 2 — Xi = s 2 e s (xi). This integral 
indeed gives the symplectic area 



a) 



AS = U2e M (xi) A s 2 e s (xi) = i 2 u 2 • 



(14) 



To generalize to arbitrary Z-encounters, we imagine a 
partner orbit 7' constructed out of 7 by I — 1 succes- 
sive steps, as illustrated for a special example in Fig. 
Each step interchanges the right ports of two encounter 
stretches and contributes to the action difference an 
amount given by (|14|l . At the same time, the two pierc- 
ings points change their position as discussed in III C II 
This step-by-step process suggests a useful transforma- 
tion of coordinates. Let Sj, Uj denote the stable and 
unstable differences between the two stretches affected 
by the j-th step. Note that in contrast to §j, uj the 
index j no longer represents encounter stretches 2, . . . , I 
but steps 1, . . . , (I— 1). Now, the change of action in each 
step is simply given by SjUj. Summing over all steps, we 
obtain a total action difference 



z-i 

AS = J2 s 3 u i 



(15) 



The transformation leading from Sj , Uj to Sj , Uj is linear 
and volume-preserving. 1 Due to the elegant form of (|TTj|) . 



1 First, consider reconnections as depicted in Fig. |7Jl for 1 = 4. 
We proceed from 17k to [71 1 in I — 1 = 3 steps. In the j-th step, 
we change connections between left ports j and j + 1, and right 
ports 1 and j + 1. Recall that stable and unstable coordinates 
of piercing points are determined by the left and right ports, re- 
spectively. Thus, the separation between the stretches affected 
has a stable component Sj = Sj+i — Sj and an unstable com- 
ponent Uj = ttj+i — iti. The Jacobian of the transformation 
s, u — > s,u is equal to 1. All other permissible reconnections can 
be brought to a form similar to Fig. |7Jl (albeit with different 
I), by appropriately changing the numbering of stretches; hence 
they allow for the same step-by-step procedure. 



1«- 
29- 

4«- 
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FIG. 7: Steps from connections in 7 (depicted in a) to those 
in 7' (shown in d), in each step interchanging right ports of 
two encounter stretches. 



it will be convenient to use s,-, Uj rather than Sj, iij in 



3 1 ■ 



defining the encounter regions, demanding all \s 
to be smaller than our bound c. Employing Q, one 
easily shows that AS remains invariant if the Poincare 
section V is shifted through the encounter. Moreover, if 
the orbits 7 and 7' differ in several encounters, the total 
action difference is additive in their contributions, and 
each is given by (|l5|l . 

At this point, we can finally appreciate that the en- 
counters relevant for spectral universality have a dura- 
tion of the order of the Ehrenfest time. The form factor 
is determined by orbit pairs with an action difference AS" 
of order H. According to our expression (|14l) for AS, the 
relevant stable and unstable coordinates are of the or- 
der \fh. The encounter duration, logarithmic in s and u, 
must consequently be of the order of the Ehrenfest time. 



D. Structures 

We want to define more precisely the notion of "struc- 
tures" of orbit pairs 7, 7'. 

(i) First of all, these structures are characterized by the 
order in which encounters are traversed in 7. We enumer- 
ate the encounter stretches of 7 in their order of traversal, 
starting from some arbitrarily chosen stretch, and assem- 
ble the labels 1 . . . L(v) in V = V(v) groups according to 
the encounters they belong to. Such a division uniquely 
defines the order in which the encounters are visited. For 
example, in an orbit pair differing in two 2-encounters the 
four stretches can be distributed among the encounters as 
(1,2)(3,4) or (1,3)(2,4) or (1,4)(2,3); each of these three 
possibilities determines a different structure. 

Some structures refer to the same orbit pair. Indeed, a 
different choice of the initial stretch in the same 7 would 
lead to a cyclic shift in the enumeration of stretches, and 
that shift may change the structure associated with 7. 
In the example of two 2-encounters, cyclic shifts may 
either leave the structure (1,2)(3,4) invariant or turn it 
into (1,4)(2,3), such that the structures (1,2)(3,4) and 
(1,4)(2,3) are physically equivalent. 

Moreover, structures are characterized (ii) by the rela- 
tive directions of the encounter stretches (i.e. =£ or 
for 2-encounters, and 3 or 2 for 3-encounters in T 
invariant systems), and (iii) by the reconnections lead- 
ing from 7 to 7'; the latter distinction is important if 
there exist several such reconnections inside the same en- 
counter set, each leading to a different connected partner. 
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E. Statistics of encounter sets 

The statistics of close self-encounters inside periodic 
orbits can be established using the ergodicity of the clas- 
sical motion. As a second ingredient, it is important to 
only consider sets of encounters whose stretches are sep- 
arated by non-vanishing loops, i.e. do not overlap. For 
example, if two stretches of different encounters overlap, 
the two encounters effectively merge, leaving one larger 
encounter with more internal stretches, see Fig. [S] The 
partners are thus seen as differing in one larger encounter, 
rather than in two smaller ones. For the more involved 
case of stretches belonging to the same encounter, see 
Appendix [Dj 

In the following, wc will consider encounter sets within 
orbit pairs (7, 7') with fixed v and fixed structure. Each 
of the V encounters of 7 is parametrized with the help of 
a Poincare section V a {a — 1, . . . ,V) crossing the orbit 
at an arbitrary phase-space point inside the encounter, 
traversed at time t a \. The orbit again pierces through 
these sections at times t a j with j — 2, . . . ,l a number- 
ing the remaining stretches of the a-th encounter. The 
first piercing may occur anywhere inside the orbit at a 
time < tn < T, T denoting the period. The remaining 
t a j follow in an order fixed by the structure at times 
tn < taj < T + tix. Each of the vi /-encounters is 
characterized by / — 1 stable and unstable coordinates 
s a j, u a j (j = — 1), which in total make for 

2 J2i >2 {l — = 2(£ — V) components. If V a is shifted 
through the encounter, the stable and unstable coordi- 
nates change while the contributions to the action differ- 
ence AS a j = s a jU a j remain invariant. 

We proceed to derive a density Wt(s,u) of phase- 
space separations s, u. To understand the normal- 
ization of wt(s,u), assume that it is multiplied with 
Ylaj 3(AS a j — s a jU a j) and integrated over all s, u. The 
result will be the average density of partners 7' per one 
orbit 7 such that the pair (7,7') has the given structure 
and action difference components AS a j. Averaging will 
be carried out over the ensemble of all periodic orbits 7 
with period T in a given time window, assuming that the 
contribution of each orbit is weighted with the square of 
its stability amplitude. According to the equidistribu- 
tion theorem this ensemble is ergodic yielding the same 
averages as integrating over the energy shell with the Li- 
ouville measure. 

We need to count the piercings through Poincare sec- 
tions parametrized by stable and unstable coordinates. 
Due to ergodicity, the expected number of such piercings 
through a given section in the time interval (t, t + dt) with 
stable and unstable components in (s,s + ds) x (u, u + du) 
is equal to ds du dt/Q, i.e. corresponds to the uniform 
Liouvillc density l/f2. 

In fact, we need the number of sets of L — V pierc- 
ings through our sections V a occurring in time inter- 



vals {taj : t/Q 



dt c 



J 



,l a , with stable and un- 



stable coordinates inside (s a j,s a j + ds a j), [u 




FIG. 8: Two 2-encounters (marked by boxes) overlap in one 
stretch and thus merge to a 3-encounter (solid bold arrows) 



by p T {s,u) d 



]L-V 



t. The uniform Liouvillc 



density carries over to the coordinates s a j , u a j since the 
transformation from s,u to s,u is volume-preserving; so 
we may expect pr(s,u) equal to l/£l L ~ v . 

However, recall that we are only interested in encoun- 
ters separated by non- vanishing loops. To implement 
that restriction, we employ a suitable characteristic func- 
tion 0t(s,u,£) which vanishes if the piercings described 
by s, u and t correspond to overlapping stretches, and 
otherwise equals 1. We thus obtain 



pr(s,u,t) = &T(s,u,t) 



1 



L-V 



(16) 



(Actually the duration of the connecting loops must 
be not just positive but also larger than all classical re- 
laxation times describing correlation decay, to guarantee 
the statistical independence of the piercings. However, 
that classical minimal loop length t c \ is not worth fur- 
ther mention since it is vanishingly small compared to 
the Ehrenfest time, the smallest time scale of semiclassi- 
cal relevance.) 

Proceeding towards wt(s,u) we integrate over the 
L — V piercing times t a j, j > 2, still for fixed Poincare 
sections V a - The integral yields a density of the stable 
and unstable components s, u of L — V piercings, reck- 
oned from the V reference piercings. To finally get to wt, 
we must keep track of all encounters along the orbits in 
question. To that end we have to consider all possible po- 
sitions of Poincare sections and thus integrate over the 
times t a i (of the reference piercings) as well. Doing so, 
we weigh each encounter with its duration t" nc , since we 
may move each Poincare section to any position inside 
the duration of the encounter. In order to count each en- 
counter set exactly once, we divide out the factors t" nc , 
and thus arrive at the desired density 



w T (s,u) = 



f d L tQ T (s,u,t) 



n 



fa ClL-V 
a L cnc J ' 



(17) 



du a j), j 



1; that number will be denoted 



It remains to evaluate the L-fold time integral in i|17[l . 
The integration over tu runs from to T; it will be done 
as the last integral and then give a factor T. The L — 1 
other t a j must lie inside the interval (iii,tii + T) and 
respect the ordering dictated by the encounter structure. 
Moreover, subsequent encounter stretches must not over- 
lap. Thus, the time between the piercings of two subse- 
quent stretches must be so large as to contain both the 
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head of the first stretch and, after a non-vanishing loop, 
the tail of the second stretch. (Given T-invariance, we 
rather need to include the tail of the first stretch, if it 
is time-reversed w. r. t. the earliest stretch of its en- 
counter; likewise the second stretch may also participate 
with its head.) These minimal distances sum up to the 
total duration of all encounter stretches ^ Q l a t" nc , since 
each stretch appears in this sum once with head and tail. 

The minimal distances effectively reduce the integra- 
tion range, as we may proceed to a new set of times t a j 
obtained by subtracting from t a j both t\\ and the sum of 
minimal distances between t\\ and t a j. The t a j just have 
to obey the ordering in question, and lie in an interval 
(0, T — J2 a lat" nc ), where the subtrahend is the total sum 
of minimal distances. We are thus left with a trivial inte- 
gral over a constant. Perhaps surprisingly, the resulting 
density 



are very small. 2 Since only the action difference is dis- 
criminated by the small quantum unit we may simplify 
the double sum (JSJ as 

K{t) = -L (y, \Ar\Wi-W6 (rT H - T 7 )\ . 

\7,7' / 

(20) 

The summation over 7 is evaluated using the rule of Han- 
nay and Ozorio de Almeida (|10|l . The diagonal pairs con- 
tribute kt, with k = 1 in the unitary and k = 2 in the 
orthogonal case. The sum over partners 7' differing from 
7 in encounters can be performed with the help of the 
density P$(AS), 



K{t) = kt + ktJ2 / dA $ Pv(AS)e l 



iAS/h 



WT[S > U) (L-iy.UA^-y 



(18) 



depends only on v but not on the structure considered, 
and that fact strongly simplifies our treatment. 

The number P$(AS)dAS of orbit pairs with given v 
and action difference within (AS, AS + dAS) now reads 

Pt(AS)dAS = dAS^l x 



(21) 



The factor n in the second member is inserted since apart 
from the partner orbits considered so far, time-reversal 
invariance demands to also take into account their time- 
reversed versions; the factor r comes from the sum rule, 
setting t = Substituting ljT§|) for P$(AS) we get 

K(t) = kt + x 



d L - y sd L - v u5 AS 



aj 



where N(v) is the number of structures existing for the 
given v. Multiplication by N(v) is equivalent to summa- 
tion over all structures belonging to the same v, since Wt 
is the same for all such structures. The denominator L 
prevents an overcounting. To understand this, remember 
that one encounter stretch was arbitrarily singled out as 
"the first" and assigned the piercing time t\\. Each of the 
L possible such choices leads to a different parametriza- 
tion by s, u of the same encounter set, and may also lead 
to a different structure. The integral over s, u in lTH?|) 
includes the contributions of all equivalent parametriza- 
tions, and this is why the factor L must be divided out. 



F. Contribution of each structure 

To determine the spectral form factor, we have to eval- 
uate the double sum over periodic orbits 7, 7' in JSJ. In 
doing so, we will account for all families of orbit pairs 
whose members are composed of loops similar up to time 
reversal, i.e. both "diagonal" pairs and orbit pairs dif- 
fering in encounters. We assume that these are the only 
orbit pairs to give rise to a systematic contribution (an 
assumption that will be further discussed in the conclu- 
sions). For the pairs 7,7' related to encounters not only 
the action difference S 7 — Sy but also the difference of 
the stability amplitudes and the difference of the periods 



d L - v sd L - v u WTM e*S^W. (22) 
L 

SajUaj )wt(s,u) (19) Here, the orbit pairs with fixed v, structures, and sepa- 
rations s,u appear with the weight N(v) WT ^ W ^ . 

The integral over s and u, multiplied with kt, yields 
the contribution to the form factor from each struc- 
ture associated to v. The integral is surprisingly sim- 
ple to do. Consider the multinomial expansion of (T — 
J2 a ^a^cnc) L1 m our expression itTB|) for the density 
Wt(s,u). We shall show that only a single term of that 
expansion contributes, the one which involves a product 
of all t" nc and therefore cancels with the denominator, 



,.contr 7^ 



(L-l)l. rph-V-l 

(L-v-iy. 



n«(-u 



L 



h(v) 



T 



li n L - v 

L-V 



(-i) y n^ 

L(L-V-l)\ 



(23) 



Due to the cancellation of i" nc , w^ ntI does not depend 
on the stable and unstable coordinates and therefore the 



2 As shown in list , the quantities determining Aj (the period, the 
Lyapunov exponent and the Maslov index of 7) may be written 
as integr als of time- reversal invariant quantities along the orbit; 
see also ll-H ll (1 l28l for the Maslov index. Since 7' locally al- 
most coincides with 7 up to time reversal, we have Ay ~ Ay . 
For the case of the Hadamard-Gutz willer m odel, a more careful 
treatment of these points is given in l2Sl l.'itl . 
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remaining integral over s and u is easily calculated, 

„.,contr 



KTh{v) 



L-V 



d L ~ v s d L - v ue^^ s ^ u " 



-> k/i(u) r L - y+1 ; (24) 
we have just met with the (L— V)-th power of the integral 



The n-th term in the series if(r) = kt + Yl n >2 ^n r ™ i s 
exclusively determined by structures with v{v) = L(v) — 
V(v) + 1 = n. It will be convenient to represent K n as 3 



K„ 



(n-2)! 
iV(iT) = iV(w) 



^(-i) v n^ ! 



(29) 
(30) 



ds e isu/fi -> 2irh 



(25) 



and used 27rfi,£ = t5- = t. In the semiclassical limit, 
the contributions of all other terms in the multinomial 
expansion vanish for one of two possible reasons: 

First, consider terms in which at least one encounter 
duration i" nc in the denominator is not compensated by 
a power of t" nc in the numerator. The corresponding 
contribution to the form factor is proportional to 



n 



ds a jdu a j 6 

''enc 



(26) 



As shown in Appendix such integrals oscillate rapidly 
and effectively vanish in the semiclassical limit, as K — > 0. 

Second, there are terms with, say, p factors t" nc in the 
numerator left uncancelled. To show that such terms 
do not contribute we employ a scaling argument. Obvi- 
ously, the considered terms may only involve a smaller 
order of T than w™ ntr ; they are of order t l ~ v ~p. How- 
ever, f2 still appears in the same order ql-v ■ To study 
the scaling with H, we transform to variables s a j = 
u a j = ^=-, eliminating the ^.-dependence in the phase 
factor of l|22l) . The resulting expression is proportional 
to h L ~ v due to the Jacobian of the foregoing transforma- 
tion, and proportional to (In h) p due to the p remaining 
encounter durations ~ In h. Together with the factor r 
originating from the sum rule, the corresponding contri- 
bution to the form factor scales like 



tT 



x 



L-V-p 

In ft 



L-V 



(lnft)P cx t 



rpL — V — p 



T 



L-V 
H 



L-V+l-f 



(27) 



and thus disappears as h — > 0, Tjj oc hr l — > oo. 

Therefore, the contribution to the form factor arising 
from each structure with the same v is indeed determined 
by l|24p. Remarkably, this result is due to a subleading 
term in the multinomial expansion of wt(s,u), originat- 
ing only from the small corrections due to the ban of 
encounter overlap. 

The calculation of the form factor is now reduced to the 
purely combinatorial task of determining the numbers 
N(v) of structures and evaluating the sum 



K(t) 



K2_^N{v)h(v) 



-L-V+l 



(28) 



III. COMBINATORICS 
A. Unitary case 

1. Structures and permutations 

To determine the combinatorial numbers N(v), first for 
systems without T-invariance, we must relate structures 
of orbit pairs to permutations. 

Most importantly, we require the notion of cycles |3l| . 
We may denote a permutation of I objects (say the nat- 
ural numbers 1,2,.. J) by {a — > P(a),a = 1 . . .1} or 

P = ( P ^x)p(2) p\i)j- Ari alternative bookkeeping 
starts with some object ai and notes the sequence of 
successors, a\ — > 0,2 = P(a>i) — > 0,3 — P{a>2) ■ ■ ■', if that 
sequence first returns to the starting object a\ after pre- 
cisely I steps one says that the permutation in question 
is a single cycle, denotable simply as (ai, 02, . . . ai). A 
cycle is defined up to cyclic permutations of its member 
objects. The number I of objects in a cycle is called the 
length of that cycle. Obviously, not every permutation is 
a cycle. A more general permutation can be decomposed 
into several cycles. 

We now turn to applying the notion of cycles to self- 
encounters of a long periodic orbit 7 and its partner 
orbit (s). We first focus on an orbit pair differing in 
a single Z-encounter. This encounter involves I orbit 
stretches, whose entrance and exit ports will be labelled 
by 1,2, ... ,1. Inside 7 the j-th encounter stretch con- 
nects the j-th entrance and the j-th exit; the permuta- 
tion defining which entrance port is connected to which 

exit port thus trivially reads = ^'"jj. A partner 

orbit 7' differing from 7 in the said encounter has the 
ports differently connected: The j-th encounter stretch 
connects the j-th entrance with a different exit P e nc(j)- 
This reconnection can be expressed in terms of a different 

permutation P cnc = ( Pcn ^ (1) Pc[ f (2 ) '.['. pJ c (i)J> e S- recon- 
nections as in Fig. 0^, are described by the permutation 
-Pcnc = (264513)' N°t e that we refrain from indexing 
the latter permutation by a superscript 7'. 



3 We slightly depart from the notation in L'(l . where N(v) was 
defined to include the denominator (n — 2)!. 
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A permutation P en c accounting for a single ^-encounter 
is a single cycle of length I, e.g. (1,2,6,3,4,5) in the 
above example. If it were multiple-cycle, the encounter 
would effectively fall into several disjoint encounters. For 
example, Fig. 0}d visualizes a permutation with three 
cycles (1,2), (3,4), and (5,6). As already mentioned, 
reconnections only take place between stretches 1 and 2, 
stretches 3 and 4, and stretches 5 and 6, which thus have 
to be considered as three independent encounters. 

If 7 and 7' differ in several encounters, the connections 
between entrance and exit ports are reshuffled separately 
within these encounters. The corresponding permutation 
P enc then has precisely one /-cycle corresponding to each 
of the vi Z-encounters, for all I > 2, the total number of 
permuted objects being L = J2i>2 ^ Vl - 

We also have to account for the orbit loops. The a-th 
loop connects the exit of the (a — l)-st encounter stretch 
with the entrance of the a-th one. These connections can 
be associated with the permutation Pi 00 p = ^231) = 

(1,2, ... ,L) which obviously is single-cycle. The order in 
which entrance ports (and thus loops) are traversed in 
7 is then given by the product P 7 = PioopPdc = -Pioop- 
This product is single-cycle - as it should be, because 7 
is a periodic orbit and hence returns to the first entrance 
port only after traversing all others. 

Similarly, the sequence of entrance ports (or, equiva- 
lently, loops) traversed by 7' is represented by 

P = PloopPmc- (31) 

with the same Poop as above. We must demand P to be 
single-cycle for 7' not to decompose to a pseudo-orbit. 

We shall denote by M. (v) the set of permutations P enc 
(representing intra-encounter connections) which have vi 
Z-cycles, for each I > 2, and upon multiplication with 
Poop yield single-cycle permutations l|31|l . These per- 
mutations Penc are in one-to-one correspondence to the 
structures of orbit pairs defined in III Dl i.e. determine 
how the encounter stretches are ordered, and how they 
are reconnected to form a partner orbit. The number of 
elements of A4(v) is thus precisely the number N(v) of 
structures related to v. 



2. Examples 

The numbers N(v) can be determined numerically, by 
generating all possible permutations P en c with v\ Z-cycles 
and counting only those for which P = PioopPmc is 
single-cycle. The P en c's contributing to the orders n = 3 
and 5 of the spectral form factor are shown in Table HJ 

Interestingly, no qualifying P on c's exist for even L — 
V — 1 = n. For example, the only candidate for n = 2 
would be P C nc = (21), describing reconnections inside 
an encounter of two parallel orbit stretches =C . However, 
the corresponding partner decomposes into two separate 
periodic orbits (corresponding to the cycles (1) and (2) 



order 


V 


L 


V 


N(v) 


N(v) 


contribution 


r 3 


(2) 2 


4 


2 


1 


1 


It 3 




(3) 1 


3 


1 


1 


-1 


-It 3 















Or 3 


r 5 


(2) 4 


8 


4 


21 


42 


7r 5 




(2)W 


7 


3 


49 


-84 


-14r 5 




(2)W 


6 


2 


24 


32 


16 5 

3 r 




(3) 2 


6 


2 


12 


18 


3r 5 




(5) 1 


5 


1 


8 


-8 

















Or 5 



TABLE I: Permutations, and thus families of orbit pairs, giv- 
ing rise to orders r 3 and r 5 of the form factor, for systems 
without T-invariance. We represent v by (2)" 2 (3)" 3 .... The 
order of each contribution is given by n = L — V + 1. We see 
that contributions add up to zero for odd n, whereas there 
are no permutations for even n. 




FIG. 9: Reconnections inside a parallel 2-encounter yield a 
pseudo-orbit decomposing into two separate periodic orbits. 

of P = Poop Pone = (12))' see Fig- ED The same happens 
for all other permutations with n even. 4 

For n odd, the individual numbers N(v) and N(v) do 
not vanish. However, we see in Table [Q that the N(v) 
related to the same n sum up to zero. That remarkable 
cancellation, a non-trivial property of the permutation 
group to be discussed below, is the reason why all off- 
diagonal contributions to the spectral form factor vanish 
in the unitary case. For example, the r 3 term is deter- 
mined by P C nc = (3412) describing reconnections in- 
side two 2-encounters of parallel stretches =t (case ppi 
in 19]), and P cn c = (231) describing reconnections in- 
side a 3-encounter of parallel orbit stretches ^ (case 
pc in 19]). The respective contributions, t 3 and — r 3 , 
mutually cancel; Fig. displays the orbits. 



4 The proof is based on the parities of the permutations involved. 
A permutation is said to have parity 1 if it can be written as 
a product of an even number of transpositions, and to be of 
parity —1 if it is a product of an odd number of transpositions. 
Parity is given by (— 1) L ~ v , where L is the number of permuted 
elements and V the number of cycles, and the parity of a product 
of permutations equals the product of parities of the factors. 
Since P and Pi oop both consist of one single cycle, they are of the 
same parity. Therefore, P = PioopPcnc implies that all allowed 
Penc need to have parity 1, i.e. n = L — V + 1 must be odd. 
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1,2, 



ppi 



FIG. 10: Two families of pairs of orbits differing in paral- 
lel encounters; both exist for systems with and without T- 
invariance; each contributes to r 3 , but the contributions mu- 
tually cancel. For labels see text. 



3. Recursion relation for N(v) 

We now derive a recursion formula for N{v), imagin- 
ing one loop (e.g. the one with index L) of an orbit re- 
moved and studying the consequences on the encounters. 
We shall reason with permutations but the translation 
rule cycle — > encounter yields an interpretation for or- 
bits. Readers wanting to skip the reasoning may jump 
to the result 

As a preparation, let us introduce a subset M(v, I) of 
M(v) such that the largest of the permuted numbers, 
i.e. L(v) = J2 k kvk belongs to a cycle of length / (it is 
assumed that vi > 0). The full set can be obtained by 
applying to this subset all L(v) possible cyclic permu- 
tations. In fact, we thus get the set M(v) in Ivi copies, 
since cyclic permutations shifting the element L(v) inside 
an /-cycle or between different /-cycles with the same I 
leave the subset A4(v,l) unchanged. Consequently, the 
sizes of A4(v) and A4(v,l) are related as 



N(v,l) 



Ivi 

W) 



N{v) 



(32) 



We need a mapping that leads from a given permu- 
tation Pone to a permutation Q C nc of smaller size, with 
a different cycle structure. Recall that any permutation 
P C nc £ M.(v,l) may be written as P cnc = P^pP, with 
Poop the single-cycle permutation (1,2, ...,L), and P 
single-cycle as well. Now suppose that the element L 
(corresponding to the entrance port following the L-th or- 
bit loop) is deleted from the cycle representations of both 
P 00 p and P. We thus obtain two new single-cycle permu- 
tations Qioop and Q, acting on the numbers 1, 2 . . . , L— 1. 
Here, Qioop is simply given by Qi oop = (1,2, 1), 
and Q differs from P only by mapping the predecessor 
of L, i.e. P~ 1 (L), to the successor of L, i.e. P(L). Let 
us now define the new "encounter" permutation Q C nc in 
analogy to (|3*T1) , 



?loop ( ' 



(33) 



, L — 1 as 



Qenc(a) 



P cnc (i) ifa = P e " 1 c (L-l) 
L-\ ifa = P" 1 c (L) 
Psnc(a) otherwise . 



(34) 



To verify this, recall that Qi 00 p differs from Pi oop only in 
the mapping of one number, same as for Q and P. Thus 
Qcnc acts like P on c on all but two numbers a. These 
exceptional cases, given in the first two lines of 13411 . 
are checked by carefully applying the above definitions 
of Qioop and Q as follows 

Q snc p-l(L-l) = QdLQP-iJUpiL-l) 



(*) 



Qiol P QP'HL) = Q^ op P(L) 

Plao P P (L) = Pencil 



QencP en l(^) ~ QloopQ^ ^oopC^) 



loop 



QTolpPP-'W 

; (35) 



here, we used P(L) ^ 1 for (*), since otherwise P en c 
would have a 1-cycle (i.e. P cnc {L) = P^P^L) = 

Pj~Q p l = L), and P(L) ^ L, since otherwise P would 
have a 1-cycle. To check (**), we need P _1 (l) 7^ L (since 
P(L) ^ 1) and P- X {1) ^ P-\L). 

We need to connect the cycle structures of Qcnc and 
Ponc- Let us first consider the case that the element 
L 1 of the permutation P onc belongs to a different 
cycle than L, say a fc-cycle. Hence, P enc has the form 

Pcnc = [. • -](L - I, a,2, a 3 , . . . a,k)(L, 6 2 , 63, ■ ■ ■ h) (36) 

where the two aforementioned cycles are written in round 
brackets, and [. . .] stands for all other cycles. Then 
Qcnc differs from P onc by mapping P e ~ c (L - 1) = a k to 
P en c(£0 = b 2 , and P e ~ c (L) = bi to L — 1. It follows that 
the k- and /-cycles of P on c merge to a (k + 1 — l)-cycle of 



Qc 



[. ..](£— 1, a 2 , as,... a k ,b 2 , b 3 



(37) 



where [. . .] is the same as in (|36H . Compared to P en c, Qenc 

has one /c-cycle and one /-cycle less, but one additional 

(k + l— l)-cycle. The changed cycle structure with v k — > 

v k - l,uj — >• vi - l,Vk+i-i — > Vk+1-1 + 1 will be denoted 
as $[k,l->k+i-i] 4 In genera^ ^ai,...,a m -/3 1 ,...,ft l ] ) m > q, 

n > denotes the vector obtained from v if we decrease 
all v ai by one, increase all Vp { by one, and leave all other 
components unchanged; if no /3, appear on the r. h. s., 
no components of v are increased. 

The permutation Q C nc thus belongs to the subset 
A4('ifi k ' l ^ k+l ~ 1 ^ , k+l—1) since the largest permuted num- 
ber L — 1 belongs to a cycle with the length k + I — 1. 
Each P enc with the structure l|3T)|) (k and / fixed) gener- 
ates one member of this subset. Vice versa, for fixed k 
the Q cnc given in l|37|) uniquely determines one P enc as 
given in l|36p . Hence, there are 



The Qenc thus obtained acts on the elements a 



[k,l— >k+l-l 



\k + l-l) 



(38) 
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members of M(v,l) structured like l|36|l . Physically, the 
present scenario corresponds to the merger of a fc- and 
an /-encounter into a (fc + I — l)-encounter, by shrinking 
away an intervening loop. 

We now turn to the second scenario where L and L — 1 
belong to the same /-cycle of P cnc - If L follows L — l 
after m iterations (i.e. L — P™ C (Z/ — 1), 1 < m < I — 2) 5 , 
the permutation P cnc is of the form 

-Pcnc = [...](L— l,a,2, 0,3, . . . a m ,L, a m+2 , ...,ai). (39) 

According to l|34|) . Q e nc differs from P cnc by mapping 
) P enc = a m+2 and mapping P e „ c (£) = 



ii^uiumj nnsy ■ ^tfenc 

P^nl( L - 1) = a; to Penc = a„ 

a m to L — 1; Qcnc thus reads 



Qo 



[. ..](£- 1, a 2 , a 3 , . . . a m )(a m+2 , . . . , a/) ; (40) 



the /-cycle of P cnc is broken up into 2 cycles, with 
the lengths m and I — m — 1. Since the largest num- 
ber L — 1 is included in an m-cycle, Qcnc belongs to 
M(^ m ' l - m - 1 \m). 

In contrast to the first scenario, there are typically 
several P enc producing the same Q e nc- Indeed Q40[) would 
not only result from i|39|) . but also from all I — m — 1 
permutations P cn c obtained by cyclic permutation of the 
last elements a m+ 2, in (|39|l . Besides, [. . .] in P enc 

contains w;_ m _i cycles of length l—m—1. If we transpose 
the content of one of these cycles with the subsequence 
a m+2, ■ ■ ■ ,cii in (|39|) . the resulting P cnc will lead to the 
same Qcnc- Thus, for each m, the subset of elements 
Pcnc S M(v, Z) structured like (|3T?|l is (I — m— l)(u/_ m _i + 
1) times larger than M.{v [ - l ^ m,l ~' m ~ 1 \m) 1 i.e. it has the 
size 



(i-m-l)(u,_ m _i + l)iV(#- 



.l — m— 1 



l,m). (41) 



We have now decomposed M,(v, I) into several subsets 
of size 1] , k + I - 1), k > 2, and further sub- 

sets of size (l-m - l)(vj_ m _i + l)^^" 1 ''-™" 1 ], m), 
m = 1, . . . , / — 1. The size of A^u, Z) thus reads 

N(v, I) = V jv^fc.'-fc+i-i] ) + j _ i) 



fc>2 



(42) 



i-2 



(l - m - l)(«i_ m _i + lJJV^™' 1 "™- 1 ]^) . 



We may rewrite the latter equation using N(v, I) — 



l -fN(v) and N(v) = N(v) y 1; L( ^'' ' , to get 

fc>2 



(-i) v U,i v ' 



l<m<l-2 

0; 



(43) 



5 Note m = I — 1 is excluded: otherwise P would have a 1-cycle due 
to L = P^iL-l) = P-I(L-I) = P-^ioopCL-l) = P- 1 ^)- 



Eqs. i|42l43|) are the general recursion relations in search. 

[k,i-»k+I-l] 



(Note that Vk+i- 



1 



v 



k+l-l 



Of course, the fcth 



summand vanishes if there are no fc-cycles present, i.e. if 
Vk = and thus formally v ^' l ~* k+l x l = —I s ). 

To determine the form factor for systems without time- 
reversal invariance, we only need the special case I = 2. 
In this case, our recursion strongly simplifies, 



v 2 N(v) 



k>2 



) = 0, (44) 



since only the first of the two above scenarios is possible. 
That is, a 2-cycle may only merge with a fc-cycle to form 
a (fc + l)-cycle, but not split into two separate cycles. 
Recall that {tf fe . 2 — ►fc+i] i s obtained from v by decreasing 
both Vk and «2 by one, and increasing Vk+i by one. 



4- Spectral form factor 

We had expressed the Taylor coefficients of the form 
factor as a sum over the combinatorial numbers N(v), 



Kn 



1 



u{v)—n 



(n-2)! 



N(v), n>2, 



(45) 



sec <|29|) . where the sum runs over all v with vi = which 
fulfill v{v) = L(v) — V(v) + 1 = n. Our recursion relation 
for N(v) now translates into one for K n , albeit a trivial 
one in the unitary case, implying that all K n except K\ 
vanish. (Alternatively, one may use a rather involved 
explicit formula for N(v) |32^-) 

To show this, consider the recursion l|44|) for N(v) and 
sum over v as above 



E 



Vn 



a> ( ,7) + Y 4 fc ; 2 r fe+11 fcA(^ fe ' 2 ^ fe+11 ) 1=0. 

k>2 



(46) 

Each of the sums E^ 71 v[ k ^ k+l] AT(#.2^fe+i] j may 

be transformed into a sum over the argument of N, i.e. 
t7' = ^[fc,2^/£+i] _ Y) ne jy^^f) = n ^ we a i so nave ^(v') = n, 

since going from v to decreases both L and V by one 
and thus leaves v invariant. Given that by construction 
we must have v' k+1 > 1, the sum over v 1 extends over 
all v' with v(v') — n and v' k+1 > 1. However, the latter 
restriction may be dropped because due to the prefactor 
v k +i t crms with v' k+1 = do not contribute. Conse- 
quently, the sum may be simplified as 



y{v)—n 

E < 



y(v )— n 



E "k+i^)- 



(47) 

Applying this rule to all terms in H46fl and dropping the 
primes, we obtain 



v{v)—n 



E (t'2+E«fe+lfc)^V( 



v) =0. 



(48) 



fc>2 



13 



Since the term in brackets is just J2i>2 v lQ~^) — L—V — 
n — 1 we have 

u(v)—n 

(n-1) J2 N(v) = (n-l)\K n = 0, n>2. (49) 



7 



7 



entrance 1 



exit 1 



exit 1 



entrance 1 




exit 2 entrance 2 




entrance 2 exit 2 



We see that all Taylor coefficients except K\ vanish: 
orbit pairs differing in encounters yield no net contribu- 
tion; the diagonal approximation exhausts the small-time 
form factor in full. GUE behavior is thus ascertained. 



B. Orthogonal case 

1. Structures and permutations 

In T-invariant systems the partners of an orbit 7 may 
involve loops of both 7 and its time-reversed 7. To cap- 
ture all partners 7' of 7 in terms of permutations the 
permuted objects must refer to both 7 and 7 and thus be 
doubled in number compared to the unitary case. Each 
permutation will describe simultaneously both 7' and 7'. 

We number the entrance and exit ports of self- 
encounters of 7 in their order of traversal, 1, 2, . . . , L, 
such that the a-th encounter stretch leads from the a-th 
entrance to the a-th exit port; see Fig. II II for the exam- 
ple of a 2-encounter. The time-reversed orbit 7 passes 
the same ports as 7, but with opposite sense and en- 
trance and exit swapped. The ports of 7 are labelled by 
1,2, ... ,L, such that the exit port a of 7 is the time- 
reversed of the entrance port a of 7, and entrance port 
a of 7 is the time- reversed of the exit port a of 7, again 
compare Fig. 1111 Consequently, inside 7 stretch a leads 
from entrance a to exit a. 

The intra-encounter connections of 7 and 7 are 
represented by the trivial permutation P^ nc = 

(12 iT3 l) which maps each entrance (upper line) 
to the following exit (lower line). The loops are associ- 

/it) I"OQ 1 \ n 

ated with Pi 00 p = I 23 "' 1 T 2 l ] ' smce if 0110 1°°P °f 7 
leads from the exit of the (a— l)-st stretch to the entrance 
of the a-th one, its time-reversed must go from exit a to 
entrance a — 1. Finally, P 7 — Poo-pP^c — P\ 00 p sped- 



r 

entrance 1 entrance 1 



T 

exit 1 exit 1 



fics the ordering of entrance ports along 7 and 7. That 
P 7 has two cycles (1,2,..., L) and (L, L — 1, . . . , 1), one 
each for 7 and 7. 

The reconnections leading to 7', 7' are described by a 
permutation P enc determining the exit port connected to 
each entrance. In the example of a Sieber/Richter pair, 
Fig. ^] the partner 7' connects the entrance 1 of 7 to 
the exit 2 of 7, and the entrance 1 of 7 to the exit 2 of 7. 
Including the connections in the time-reversed partner 

7', we write P cnc = ^Y2i)- Note that the sequence 
of columns in P onc may be ordered arbitrarily. We shall 
mostly order such that the first lines in P enc and P^ nc 
coincide; columns describing 7' and 7' may thus become 
mutually interspersed. 




exit 2 



exit 2 




entrance 2 



entrance 2 



FIG. 11: Entrance-to-exit connections for a Sieber/Richter 
pair, in an orbit 7, its time-reversed 7, and the partners 7', 7'. 



T-invariance imposes a restriction on P cnc : If a stretch 
connects entrance a to exit b, the time-reversed stretch 
must connect the entrance b of 7 to the exit a. Thus, 
if P enc maps a to 6, it has to map b to a, with a, b 
standing for elements out of 1, . . . , L, 1, . . . , L (we define 
a = a). This restriction on P enc will be referred to as 
"T-covariance" . 

It follows that if (a%, a%, . . . , oj_i, a{) is a cycle of 
P onc , then so is its "time-reversed" , (a7, aj_i, . . . ,02, aT). 
These two cycles may not be identical. Indeed a cycle 
coinciding with its time-reversed would have the form 
(ai, . . . , afe,Ofc, . . . ,ai) , k — 1/2] such cycles are not al- 
lowed, since the entrance port Ofc and the exit port a*, 
coincide in configuration space and thus may not be con- 
nected by an encounter stretch. 

We see that given time reversal there must be a pair 
of twin ^-cycles of P onc associated with each Z-encounter; 
an encounter associated with more than one pair of cy- 
cles would decompose into several independent ones. In 
general, each cycle in a pair describes stretches both of 
7' and of 7'; only in case of all stretches nearly parallel 
one cycle refers exclusively to 7' and the other to 7'. 

The final restriction on the permutations P cn c is anal- 
ogous to the one encountered in the unitary case. To ob- 
tain two connected partner orbits 7' and 7', we now have 
to demand the permutation P = PioopPmc to consist of 
only two L-cycles, listing the entrance ports in 7' and 7', 
respectively. (The second cycle can also be interpreted 
as the list of exit ports of 7', time reversed and written 
in reverse order. Since an entrance a t is connected by a 
loop to the exit b t = P^ p (ai), the two cycles of P read 

(ai,a 2 , . . . ,a L ) and {b L ,b L -i, ■ ■ -M)-) 

To summarize, we consider the set M(v) of permuta- 
tions P Gnc acting on 1, 2, . . . , L, 1, 2, . . . , L which (i) are 
time-reversal covariant, (ii) have vi pairs of Z-cycles for 
all I > 2, and (iii) lead to a permutation P = Pi oop Penc 
consisting of two cycles as above. Each element P cnc of 
the set Ai(v) thus stands for one of the "structures" in- 
troduced in III Dl We need to calculate the number N(v) 
of elements of M. (v) . 
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order 


V 


L 


V 


N(v) 


N(v) 


contribution 


r 2 


(2) 1 


2 


1 


1 


-1 


-2r 2 












-1 


-2t 2 


r 3 


(2) 2 


4 


2 


5 


5 


10r 3 




(3) 1 


3 


1 


4 


-4 


-8r 3 












1 


2r 3 


4 

T 


(2) S 


6 


3 


41 


164 

:-i 


164 4 
3 




(2) 1 (3) 1 


5 


2 


60 


72 


72r 4 




(4) 1 


4 


1 


20 


-20 


-20r 4 












8 




T 5 


(2) 4 


8 


4 


509 


1018 


1018 5 
.3 T 




(2) 2 (3) 1 


7 


3 


1092 


-1872 


-624r 5 




(2) 1 (4) 1 


6 


2 


504 


672 


224r 5 




(3) 2 


6 


2 


228 


342 


114r 5 




(5) 1 


5 


1 


148 


-148 


148 5 
3 












12 


4r 5 



TABLE II: Permutations, and thus families of orbit pairs, 
giving rise to orders r 2 to r 5 of the form factor, for systems 
with T-invariance; notation as in Table|l] The results coincide 
with the predictions of RMT for the GOE. 



2. Examples 




api 

FIG. 12: Pairs of orbits existing only for systems with T- 
invariance and giving all of r 3 . For labels see text. Two 
further families do not require T-invariance, see Fig. 1101 



family api), or (iii) two antiparallel ones (P cnc = 

^4 3 214 3 2 l) ana r <=nc - ^2 14 3 2 1 4 3^' Iamll y aas 

in 0]). Again, equivalent permutations differ by cyclic 
permutations (of now 1,2,3,4 as well as 1,2,3,4), i.e. in 
which of the stretches is assigned the number 1. 



Again, the numbers N(v) can be determined by numer- 
ically counting permutations. From the results shown in 
Table ILTI we see that indeed the form factor of the Gaus- 
sian Orthogonal Ensemble is reproduced semiclassically. 

The r 2 contribution comes from pairs of orbits differ- 
ing in one antiparallel 2-encounter Q- We have already 
shown that the corresponding "encounter permutation" 

reads P onc = f^fai)- 

The t 3 contribution originates from (compare Figs. 
1101 ITS! and Ref. [19() four permutations related to 3- 
encounters and five permutations related to pairs of 2- 

encounters. The permutation P onc = (2313T2) °- e " 
scribes encounters of three parallel orbit stretches =t 
(case pc). For triple encounters in which one of the 
stretches is time-reversed with respect to the other two 
iE? (case ac), there are three related permutations, 

P cnc = ^33121!) an d its two images under cyclic per- 
mutation of 1,2,3 as well as 1,2,3; physically, the three 
latter are equivalent since they differ only in which of 
the three stretches is considered the first. (In Ref. 0] 
such equivalences were taken into account by multiplicity 
factors Afp C and Af ac -) 

Pairs of 2-encounters may either be composed of either 
(i) two parallel encounters (P cn c = (3412IIH)' ^ am " 
ily ppi) or (ii) one parallel and one antiparallel encounter 

(P - fl 2 341~ 2 3 4\ , p _ fl 2 3 4T 2 3 4\ 

v'enc ^34r234r2y cU1U ^ cnc — ^34r234t2,l' 



3. Recursion relation for N(v) 

We are now fully equipped to establish a recursion re- 
lation for N(v) in much the same way in the unitary case. 
Impatient readers may want to jump to the resulting re- 
cursion (|68() for K n . 

First of all, we recover Eq. (EH, N(v,l) = ^N(v) 
using exactly the same arguments as in the unitary case. 
From each element P enc G M. (v,l) we obtain L elements 
of M.{v) by the L possible cyclic permutations. Applying 
these cyclic permutations to the set A4(v, I) we obtain vil 
copies of M. (v) , since the vi pairs of twin ^-cycles have al- 
together vil members without overbar, and permutations 
shifting L among these members leave M(v, I) invariant. 

Choosing a permutation P cnc = P^„ p P 6 M(v, I) we 
set Qenc = QioopQ witn Qioop and Q obtained from P loop 
and P by omitting L and L — 1, and replacing L by L — 1. 
(This prescription can be interpreted as removing the 
loop leading from exit L — 1 to entrance L, and its time- 
reversed leading from exit L to entrance L — 1; hence 
the corresponding entrance ports must be removed from 
P and Poop = P 7 .) In particular Qi 00 p will have the 
form Q loop = (1,2,...,L- l)(r^T,r^2,...,T). The 
two cycles of Q fulfill the same relation as those of P 
and may thus indeed be interpreted as lists of entrance 
ports of two time reversed orbits. The resulting "en- 
counter permutation" Q cnc maps the remaining elements 
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=(a) 



1,2, ... ,L - 1,1,2, . . . ,L - 1 as 

Penc(L) ifa = P-l(L-l) 
P enc (L-l) ifa = P- 1 c (i) 
L-l_ ifa = P f gL) 

Penc(i) if = 1-1 

Penc(a) otherwise. 



(50) 



Here, the second and fourth line extend (|34[l as required 
by T invariance; the present Q C nc is indeed T covariant. 

When analyzing the cycle structure of Q cnc , we now 
have to distinguish three cases, the first two paral- 
leling the treatment of IIII A 31 Note however a fac- 
tor 2 appearing in the second case. For each Q cnc £ 
A4(if- l ~ ►"M-"*— !] ; m), there are now twice as many, 
namely 2(7 — m — l)(uj_ m _i + 1) related P cnc £ Ai(v, I) 
structured like 1|39|) . since <5 cnc also remains unaffected 
by time reversal of a m +2, ■ in (|39[1 . The second and 
fourth line in (|50|l make sure that merging or splitting of 
cycles is mirrored by the respective twins. 

A third possibility appears since the cycles involving 
L and L — 1 may be twins, and hence belong to the 
same encounter. Since the twin cycles are mutually time 
reversed there is one cycle containing both L and L — 1, 
and another one containing L and L — 1. Assume that 
inside the first cycle, the element L — 1 follows L after 
m iterations, i.e. L — 1 = P™ C (L) (with 1 < m < I — 1). 



Pcnc — [• ■ ]{L, Cl2, ■■■ , a m , L—l, a m+ 2, ■ ■ ■ , a;) 

(o7, . .. ,a m+2 , L - 1, a^T, . . . ,02, L) . (51) 

Due to Q5()jl. Qcnc differs from P onc by mapping 

Qcnc(«m+2) = &2 Qenc(a2) = &m+2 

<2cnc(a;) = ^ ~ 1 Qcnc{L - 1) = a7 (52) 

The initial pair of twin cycles of P enc is transformed to 
the following pair of twin (I — l)-cycles of 



Qenc — [• • -J(ct2, ■ • • , «m, L — 1, d/, ... , a m+ 2) 

(a m +2, ■ ■ • , au L - 1, a^T, . . . ,02) . (53) 

Given that the largest number permuted by Q C nc, i.e. 
L — 1, is included in one of these cycles, we have Q e nc £ 
7\4(tj[^' _1 l, l _ 1). Conversely, for any such Q e nc and 
each 1 < m < I — 1, there is exactly one related P cn c, 
since we may read off ai , . . . , ai from Q cnc in (|53|l and 
recombine them to form a permutation P cnc as in Q51[l. 
We thus see that each of the I — 1 subsets of M (v, I) with 
~ 1 = P™ C (L) is one-to-one to ^(-u^" 1 ], Z - 1) and 
thus has an equal number of elements. 

We have seen that A4(v, I) falls into subsets similar to 
the unitary case, time-reversal invariance making for the 
factor 2 explained above, and for I — 1 additional subsets 
of size N {0^-^,1 - 1). The various sizes combine to 



the orthogonal analogue of the recursion relation (|42() , 
N(v, I) = ^2 N(rf k ^ k+l -^ ,k + l-l) 



k>2 



1-2 



+ 53 2(Z - m - + l)N(v 

+ (l-l)N(^ l - 1 \l-l), 



— >m,l—Tn— 1] 



(54) 



which using (|32|) and the shorthand N(v) — 
N{v)- — j~ — may be written as 

viN(v) + £)(«fc+l-i + l)kN(v^ k+l -^) 



k>2 



1-2 



+ 2(«I-m-l + 



77i,l — m— 1] J\f(^y[l — — W— 1] 



-(Z-l)(^_ 1 + l)iV(u^- 1 ]) = 0; 
recall that Ufc+z-i + 1 = vf^^ +l ^. 



4. Spectral form factor 



(55) 



Similarly as in IIII A 4l we now turn the recursion rela- 
tion for N(v) into one for the Taylor coefficients K n . As 
a preparation we generalize our rule l|47|l . For all similar 
sums over v with fixed v(v) — n and Vi — we find 



f (v)— n 

£ 



y^[a 1 ,Q2,...->/3]^^-[Qi,Q 2 ,...^/3]^ 



f(if)N(if), 
(56) 

with integers > 2, /3 > 2, n' = ^(w*) = 
u (tf ai ,a 2 , . . ) = n _^2 .( a . - 1)+ (/3_ i) } an d / any func- 
tion of v' vanishing for v' j3 = 0. Eq. (|56|l follows in the 
same way as l|47|l. i.e. by switching to if = jjl ai ' a2 i— ~ *P\ 
as the new summation variable and dropping the restric- 
tion v'p > 1 (if with v'p — do not contribute due to the 
vanishing of /). One similarly shows that the foregoing 
rule holds even without any conditions on / if (3 is re- 
moved, i.e. if no new cycles are created. It is convenient 
to abbreviate the r. h. s. of (|56|1 with the help of 



v{y )—n 



(57) 



for arbitrary /; we note that K n = ( n _ 2 )i S n [1] . Thus 
equipped we turn to the three special cases I = 1,2,3 of 
our recursion relations I|54I55I) which we shall need below. 

First, the case 1=1 involves permutations with 1- 
cycles, appearing only in intermediate steps of our calcu- 
lation. If the clement L forms a 1-cycle, it may simply be 
removed from a permutation without affecting the other 
cycles, which corresponds to a transition v 



thus have N(v, 1) = N^ 1 ^) and equivalently 

Vl N(v) + L(t/ [1 - >1 )iV(ir [1 - >1 ) = . 



We 



(58) 



16 



Second, for I = 2 (and V\ — 0) the recursion (|55|) boils 
down to 



v 2 n(v) + (vt 2 r k+1]) kN(v^ k+1 ^ 



k>2 



(59) 



where only the last term is new compared to the unitary 
case. We bring it to a form free from 1-cycles by invoking 
|J5SJ and thus N(v^ 2 ^) = -L(v^)N{v^), to get 



t; 2 



i [As,2-A+l]) AJ y-/^fc,2- 1 .fc+l] 



'fc+1 



k>2 

+L{^)N{^) =0. 



(60) 



As in the unitary case, we sum over all v with v\ = and 
2/(iT) = n. The rule (|56|l and the shorthand (|57|l give 



C 2 



fc>2 



+ S n - 1 [L(v)} 



= (n-l)S n [l] + S n . 1 [L(v)} = 0. 



(61) 



A further relation is obtained by multiplying H6()(l with 
L(v) - 1 = L(#' 2 ^' ! + 1 ]) = L(z7l 2 ^l) + 1 and again sum- 
ming over v with the help of (|56|l . The resulting equation 



(n-l)5„[i(i?)]-5 n [«a] + S n _i[i(fl)(L(t/) + l)] =0 (62) 
can be simplified by l|61|) and replacing n — > n — 2, 

-(n-2)(n-3)5 n _ x [l] = 5„_ 2 b2]--5„_ 3 [L(^(L(v)+l)] . 

(63) 

Finally, we consider I = 3 (and v\ = 0) 



« 3 ^) + J2 (4 



fc>2 



) _[fe,3-^fc+2] fe ^^fe, 3 ^fc+2] 



+4A>(t/[ 3 ^ 1 - 1 ]) - 24 3 " 21 ^V(« [3 ^ 2] ) - . (64) 

The 1-cycles in the third term are eliminated using the 
identity /V^ 3 "* 1 ' 1 ]) = ±L{v^){L(v^) + l)N(v^), 
which follows by twice applying l|58|l to 

ul 3 ^ 1 . 1 ]. Sum- 
ming over v in l|64(l we find 



Sn 



V3 + ^ v k+2 k 



k>2 



+ 2S n „ 2 [L{v)(L{v) + l)} 



-2S n - 1 [v 3 ]=0. 



(65) 



This expression can be simplified using v% + 
J2k>2 v k+2k = J2k>2 v k(k - 2) = L(v) - 2V(v) = 
2(u(v) - 1) - L(v), i.e. 

2(n - 1)5»[1] - S n [L{v)] + 2S n - 2 [L(v){L(v) + 1)] 
-2S , n-i[« 2 ]=0. (66) 



Finally applying (|61() . substituting n — > n— 1, and divid- 
ing by 2, we proceed to 



71-1 



S„[l] + (n-2)5 n _i[l] 



= S n -2N - S n - 3 [L(v){L{v) + 1)] 



(67) 



Upon comparing the recursion relations (|63J) and H67J) . 
obtained for the cases Z = 2 and I — 3, we find the co- 
efficients 5 n [l] = K n and 5„_i[l] = — 2 3 ^ ! K n _i 
related as ^5„[1] = -(n - 2) 2 5„_i[l] or 



(n -!)«■„ = -2(n-2)X n _i 



(68) 



An initial condition is provided by the Sieber/Richter 
result for orbits differing in one 2-encounter, K 2 = —2. 
Thus started, our recursion yields the Taylor coefficients 



K r , 



(-2Y 



71-1 



(69) 



coinciding with the random-matrix result. Universal be- 
havior is thus ascertained for the small-time form factor 
of fully chaotic dynamics from the orthogonal symmetry 
class; the resulting series converges for r < i. 



IV. SPINNING PARTICLES AND THE 
SYMPLECTIC SYMMETRY CLASS 

We now allow for a spin with arbitrary but fixed spin 
quantum number 5. Assuming time-reversal invariancc 
we know that for integer 5 the time reversal operator T 
squares to unity, T 2 = 1, whereas for half- integer spin, 
T 2 = -1: we face the orthogonal and symplectic symme- 
try class, respectively. The semiclassical theory of spin- 
ning particles is discussed in |3^; off-diagonal terms of 
the form factor were considered in a preliminary version 
in |2lj , and for quantum graphs in [23] ■ 

The Pauli Hamiltonian reads H = Hq+I1 ■ S where S is 
the spin operator, with S 2 = h 2 S(S + 1). The vector op- 
erator B(q, p) describes the influence of the translational 
motion and external fields on the spin. The spin-orbit in- 
teraction formally behaves like 0(H) and tends to zero in 
the semiclassical limit; it is not a small quantum per- 
turbation though, since its matrix elements are infinitely 
larger than the energy spacing A ~ 0(hf), f > 1. 

The state of the system is given by a spinor with 25+1 
components and the propagator by a (25 + 1) x (25 + 1) 
matrix. In the leading order of the semiclassical approxi- 
mation the propagator consists of the scalar translational 
part which is the van Vleck propagator of the spinlcss 
system, multiplied by the spin evolution matrix d (be- 
longing to the spin-5 representation of 5/7(2)); the 
latter matrix has to be evaluated on the classical orbit (of 
the translational motion) connecting the initial and final 
point. Along such an orbit 7 the d-matrix is a function 
of the initial and final times, d = d 7 (t, to). It satisfies the 
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equation id t dy[t,t ) = B 7 (t) ■ (S/H) dy[t,t Q ) where the 
classical time-dependent vector B 7 (t) is obtained by sub- 
stituting the classical coordinates and momenta along 7 
for the operator valued arguments of B(q, p); the initial 
condition is dy(to,to) = l2S+i- Such semiclassical treat- 
ment keeps the translational motion unaffected by the 
quantum spin. The spin, however, is driven by the trans- 
lational motion. No semiclassical approximation for the 
spin itself (which would require the assumption of large 
5) is invoked. 

The full quantum nature of the spin (finite 5) notwith- 
standing, a seemingly classical manner of speaking about 
the spin is possible, due to the following fact: any matrix 
from can be parametrized by three Euler angles (e.g. 
like d{0,<j),ip) = e ^s x /h e ies x /H^s x /h^ which are time _ 
dependent for d = d(t). The angles tp(t),9(t),(j)(t) may 
also be imagined to specify the orientation of a fictitious 
rigid body in classical rotation; as done in pll I22I l33| we 
shall speak of that motion as "spin rotation" . 

We assume the translational motion chaotic as be- 
fore and require ergodicity of the combined spin rota- 
tion and translational motion. The spin rotation itself 
is then also ergodic. This means that time averages of 
the spin dependent properties over intervals longer than 
a certain relaxation time t c \ can be replaced by aver- 
ages over all d £ D^ s '; the measure to be used reads 
dfi = d<f> sin 8 d6 dip /8tt 2 . 



A. Integer spin 

The trace formula for a particle with spin j3^| gives the 
level density as a sum over periodic orbits 7 



p osc {E) ~ -i-Re V(tr d 7 )A 7 e i 
irn * — ' 



(70) 



beside the stability amplitude (including period and 
Maslov phase) and the classical action of the 7th orbit, 
the factor tr d 7 = tr dy(to + T 7 ,to) appears and reflects 
the spin evolution over a period of the translational mo- 
tion. That tr dL is independent of the initial moment t$: 
its shift leads only to a similarity transformation of d 7 . 
The form factor becomes the double sum 

K ^ = ^(E( trd 7)(trdy)A^y 



x e i(s,-s,,)/n § f rTg - T T + Ty 



(71) 



Due to the spin, the average level density and thus the 
Heisenberg time Tjj are increased by the factor 25 + 1. 
The diagonal approximation yields the sum 

^diag(r) = (]T(tr d^\Ay\ 2 S(rT H - T 7 ) ) . (72) 



Since the spin dynamics is ergodic and since we are av- 
eraging over an ensemble of orbits, the equidistribution 
theorem j2|| allows us to treat cL as random and to inte- 
grate over all matrices of the spin-5 representation 
of 5/7(2); the sum over 7 gives the usual factor T. The 
spin integral gives unity for any 5, and so we have |33| 



^diag(r) = 2r / dfi(A)(tiA) 



2r. 



(73) 



As to off-diagonal contributions from orbit pairs 7, 7', 
spin makes for two modifications compared to the pre- 
vious Sections. First, Tu contains a factor 25 + 1, such 
that the L - V factors 2nh^ in $Z%$ give (25 + 1)^ = 
(25+ l)r rather than t. Second, the contribution of each 
orbit pair comes with the factor (tr cL)(tr dy). 

To evaluate the factor (tr d 7 )(tr dy), we decompose 
the orbit 7 into L pieces by cutting it once in each en- 
counter, and represent d 7 as a product of L matrices 
describing the spin evolution over one orbit piece. In 
the orbit pairs contributing to the form factor the du- 
ration of each piece (the orbit loop + segments of the 
preceding and following encounter stretches) exceeds the 
Ehrenfest time Tg , and thus also the classical relaxation 
time t c i. Therefore, keeping in mind that we are sum- 
ming over an ensemble of orbits, we may invoke ergod- 
icity and replace the partial spin evolution matrices by 
matrices randomly chosen out of D^ S K Numbering the 
orbit pieces and the corresponding random matrices in 
the order of their traversal in 7 we may replace d 7 by a 
product AlAl-i ■ ■ - Ai, the earliest propagator matrix 
written rightmost. The orbit partner 7' consists of prac- 
tically the same pieces passed in a different order, some 
of them in the opposite direction. Hence dy can be ex- 
pressed in terms of the same L matrices At, but with 
the order suitably rearranged and with Ai — * -47 1 for 
the time reversed pieces. The expectation value of the 
trace product can now be evaluated as an integral over 
Ai, i = 1, . . . , L. Using the results of [13, one obtains 

(tr d 7 )(tr dy) 

- j(tvA L A L -i... Ai){t*Al L L AT-_\ ...Al\)f[d»(Ai) 



2S 



V25 + 1 



L-V 



(74) 



most remarkably, for an orbit pair with a given v the L- 
fold integral depends only on the difference L(v) — V(v) = 
n— 1; in particular, it is independent of the special order- 
ing of loops in the partner orbit 7' (expressed by the sub- 
scripts fci, &2, • • • > as well as of the senses of traversal 
(expressed by the exponents r\ a — ±1). 6 



If desired, these indices can be determined from the permutation 
P = Pi 00 pPonc describing the partner 7', see Section ^1111 The 
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Now the two occurrences of (25 + \) L v mutually can- 
cel and the form factor reads 

K(t) =2t + 2J2 N(v)h{v){~l) 2S{L - v 'W L - v+1 . (75) 

V 

For integer spin, (— 1) 2S = 1 whereupon we recover the 
expansion (|28|l of K(r) for the orthogonal class. 

B. Half-integer spin 

For half-integer spin, the minus sign in H75[) becomes 
relevant. Moreover, all levels become doubly degener- 
ate a la Kramers Q . With the density of levels reduced 
to half the density of states we are led to the rescaling 
K(t) — > ^K(^) [33. In this case, the form factor reads 

K ^ = 7j-X>(*W) {-l) L V+1 = i-I ln (l-r). 

V 

(76) 

To understand the final step, we compare the sums over 
v in (|76|l and (|28l) . the latter pertinent to the orthogo- 
nal class, and find K(t) = — ^Kqoe{— §)■ We have thus 
verified the random-matrix result for the Gaussian sym- 
plectic ensemble © • As predicted in [2l| the sign change 
of the argument r, which entails the logarithmic singu- 
larity of the symplectic form factor at r = 1, comes from 
the contributions of the spin integrals 174|l . 

V. RELATION TO THE u-MODEL 
A. Introduction 

The so-called sigma model [3 is a convenient 
framework for calculating averaged products of Green 
functions of random Hamiltonians. Its zero dimensional 
variant affords, in particular, the two-point correlator of 
the level density (and its Fourier transform, the spec- 
tral form factor) for the Gaussian ensembles of RMT; 
see Appendix [E] for a brief introduction. Perturbative 
implementations exist for the three Wigner/Dyson sym- 
metry classes and yield the respective spectral form fac- 
tors K(t) as power series in the time r, i.e. precisely the 
series extracted from Gutzwiller's semiclassical periodic- 
orbit theory in the preceding sections. 

The sigma model for random-matrix theory proved of 
great heuristic value for our semiclassical endeavor: We 
were led to the correct combinatorics of families of or- 
bit pairs by an analysis of the perturbation series of the 



permutation P consists of two L-cycles relating to 7' and 7' ; the 
sequence of loops fei . . . hi, in 7' is given by the appropriate cycle 
in which the elements with a bar (indicating time reversal of the 
loop) are modified like k — > (k + 1) mod L and simultaneously 
the associated superscripts r\ are set to -1. 



sigma model. The analogy of periodic-orbit expansions 
to perturbation series might prove fruitful for future ap- 
plications of periodic-orbit theory, and that possibility 
motivates the following exposition. 

Before entering technicalities it is appropriate to point 
to some qualitative analogies and differences between the 
two approaches. Very roughly, different Feynman dia- 
grams of the sigma model (both for the Wigner/Dyson 
ensembles and disorder) correspond to different fami- 
lies of (pairs of) periodic orbits, vertices to close self- 
encounters, and propagator lines to orbit loops. An im- 
portant difference lies in the point character of vertices 
and the non-zero duration, of the order of the Ehrenfest 
time Te oc Ink, of the relevant self-encounters. Of course, 
the relevant encounter durations are vanishingly small 
compared to the typical loop lengths (~ Tjj oc hr^ +1 )\ 
nevertheless, we may say that self-encounters give inter- 
nal structure to vertices. 



B. Expansions of two-point correlator and form 
factor 

The connected two-point correlator of the density of 
levels, R = (jjp — p 2 )/p 2 can be obtained from the en- 
semble averaged product of the retarded and advanced 
Green functions G±(E) = (E±i5- H)- 1 as 



, , Re 5 tr G+(E + s/2np) 5 trG-(E - s/2np) 
R ( s ) = TlPi ' 

5trG± = tr G± - trG± . (77) 

Here, the overbar denotes an average over a Gaussian 
ensemble of random matrices. The argument E (the av- 
erage energy) is suppressed. The energy difference is 
expressed in terms of the dimensionless offset s. The 
Fourier transform of R w.r.t. the offset s gives the cen- 
tral object of the present paper, the spectral form factor, 

KM = - / dsc 2isT R{s). (78) 
t J-00 

As briefly shown in an Appendix, a bosonic replica 
variant of the a-model yields the (l/s)-expansion of the 
Fourier transform of the small-time form factor as an 
integral over matrices B, 

R(s) ^Re lim r ^oi<9 s 2 s- Kr2 

2 r z 

x J d Be^^^ sl "^ BB ^\ (79) 

where on the r.h.s. the offset s must be read as s + i<5 
with 5 I 0; the matrices B, B> are r x r for the GUE and 
2r x 2r for the GOE; as before, the factor k takes the 
respective value 1 and 2 for the two classes. The essence 
of the replica "trick" is to find the foregoing integral as 
a power series in r and to isolate the coefficient of r 2 . 
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In the limit s —* oo the principal contribution to R(s) 
comes from the Gaussian factor exp[(2i/«;) trM] in the 
integrand of 1)79(1 . where M = BB^ . The remaining fac- 
tor can be expanded as 



cxp I ^jy-JtrM' 



l>2 



easily shows that the Taylor coefficients of K{r) 
K \ T \ + E«>2 K nk\ n are given by 



K„ 



L(v)-V(v) + l=n y 



(n-2)! 



E 



-V-l 



ri/>2 v i 



■N c (v) . 



(84) 



E^f?) (5>-w 



?ri/>2^ ! U 1 



v '>2 
v\ V 



V-L 



JJ(trM 



(80) 



l>2 



Here, the summation extends over integers 
V2, «3, • • • , vi, . . . each of which runs from zero to 
infinity, and we write v — (t>2, t>3, . . .) just like in our 
semiclassical treatment. The total number of traces in 
the summand v is V(v) = X)z>2 v ii anc ^ a g a i n we define 
L(v) ~ J2i>2^ Vl - The integral over B and B' in l|79|) 
may now be seen as a sum of averages like 



</(B,flt)) = J d Bf(B,B^)c 



(2i/fo)tr(BB t ) 



(81) 



We may thus write 
R(s) , 



i>2^ ! \« 



2i 



1 

t^o r 2 " s ' 
v 

„V-L 



(82) 



\Z>2 



The leading term w = corresponds to (1) = 

(const) Kr ^5 1. The respective contribution to the 

two-point correlator is — ^Re lim r ^ ;3-d 2 (s + i5)~ Kr2 = 
and thus k\t\ for the spectral form factor, 



-Re 



reproducing the diagonal part both in the unitary and 
orthogonal cases. 

For all other terms, the operations of taking the second 
derivative by s and going to the limit r -> commute, 

2 

meaning that the factor s~ KT can be disregarded. We 
shall show below that the averages of the trace products 
with non-zero v have the property 



i 



(-2i) L W 



(83) 



C. Contraction rules 

In the following, we will derive a recursion for N c (v) 
similar to the recursion in our semiclassical analysis. To 
that end, we calculate the averages of the products of 
traces in (18311 by Wick's theorem. Each average becomes, 
for the GUE, a sum of contractions of a fixed matrix B 
in one of the traces and all matrices B^; for the GOE, 
contractions with other matrices B arise as well. In all 
formulae below, X and Y must provide the traces on the 
1. h. s. with an alternating sequence of B's and B^s; then 
the same will hold for all traces on the r. h. s. . Moreover, 
the term [. . .] will stand for inert traces unchanged by the 
contraction. The GUE involves two contraction rules, 



1 



i(tr(Xy)[...]) (85a) 



(tr BX tTB^Y{. . .]) = 

(tr BXB^Y [...]) = -7y( tlX trr [...]>. (85b) 



For the GOE, two more rules arise from contractions of 
B with B (and similarly B^ with B^): 



(tr BX tvBY[. . .]) 



|<tr(xyt )M ) (85c) 



(tr BXBY[...}) = -I(tr(XYt)[...]>. (85d) 

The only possible ordering of B,B^ after contraction is 
alternation BB^ BB^ .... We may thus express all quan- 
tities in terms of M = BB^ . Each contraction reduces 
the number of M's by 1. 

The sequences X, Y in (|85b ) may be absent; then they 
must be replaced by the unit matrix 1. In particular, 
if we repeatedly invoke H85b-d) in order to reduce the 
number of M's, the final step will always be 



n 



i 



i 



(trBBt) = --((trl) 2 ) = -- K V(l) 



2i 



(86) 



in which N c (v) take positive integer values; we shall 
in fact come to interpreting N c (v) as the "num- 
ber of contractions"; the traces of M l will be 
called ^-traces, to stress the analogy with the l- 
encounters of periodic orbits. The form factor K(t) 
is now obtained by Fourier transforming. Using 
I-oo fe 2i;sr Re [r n+1 (s + i(5)-™- 1 ] = (-2) n |r| n /n! one 



Thus, in the limit 



all averages of trace products 



vanish like r or faster. Since only terms 



count, 



the contractions between the neighboring B and B^ in 
the same trace may be disregarded, unless we are dealing 
with case JUl- (Taking in ijgHbl X = 1, Y or [. . .] not 
equal to unity, the r. h. s. would be r times another 
averaged trace product and thus vanish like r 3 or faster.) 
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D. Recursion formula for the number of 
contractions 

To translate H85h,-c0 into a recursion relation for the 
numbers of contractions N c (v), let us select a trace inside 
the product in say trM z (assuming vi > 0), and 

a matrix B inside. We must contract that B with all 
other suitable matrices in the product of traces. Three 
possibilities arise, paralleling the recursion relation for 
the combinatorial numbers iii lllll 

[i] First, we take up the contractions between our se- 
lected B in trM* and all suitable matrices in some other 
trace trM fe . In the unitary case rule H85fa ,~) implies that 
one fc-trace and one Z-trace disappear while a (k + 1 — 1)- 
trace is born 



< trAf' trM k [. . .]) = _i(trM fc+I - 1 [. ..]> . (87) 

The contractions with all matrices in fc-traces trM fe 
give the same result. We thus get k(vk — Ski) contribu- 
tions like (|57Jl . where S k i is subtracted to exclude con- 
tractions between matrices inside the same trace. In 
the orthogonal case we must also invoke rule (|85b ) for 
contractions with k(vk — Ski) matrices B in traces tr_M fc 
which again all contribute identically. 

Each time, one fc-trace and one /-trace disappear and 
one (k + 1 — l)-trace is added to the trace product. The 
vector v thus changes to according to Vk — > Vk — 1, vi — > 
Vi — 1, Vk+i—i — > Vk+i-i + 1; using the same notation as 
in our semiclassical analysis we write $k,l-*k+l—\] _ The 
overall number of matrices M is decreased by 1 such 
that L — * L — 1. Each of the above contractions pro- 
vides a contribution 7V c (t/[' c ' i ^' i "+'- 1 ] here, the denom- 
inator (— 2i) in the contraction rules is compensated by 
the factor (— 2i) L ^ in the definition of the contraction 
numbers. Thus, the overall contribution to N c (v) reads 
Kk(vk~S k i)N c (^ l ~* k+l -^). 

[ii] Next, we turn to "internal" contractions between 
the selected B and all matrices B^ in the same trace 
trM', apart from those immediately preceding or follow- 
ing the selected B. As explained above, the latter con- 
tractions would increase the order in r and lead to a result 
that vanishes as r — > 0. We apply rule (JSSJd) and replace 
trM 1 by the product of two traces which together contain 
I — 1 factors M. Thus, one Z-trace disappears and two 
traces, of M m and M /_m_1 , are added, with m running 
1,2,. ..,1 — 2; the vector v then changes to #['- > "M- m - 1 ] . 
From each of these contractions, N c (v) receives a contri- 
bution N c {v^ m - l - m -^). 

[iii] For the orthogonal case rule I|85t yields further in- 
ternal contractions: Pairing the selected B with all other 
I — 1 matrices B, we obtain I — 1 times trM' -1 . Each 
time v is thus replaced by v^ l ~* l ~^. Altogether, we gain 
the contribution (I - 1)N C (^ 1 -^). 

Summing up all contributions, we arrive at the recur- 



rence for N c , for any I with vi > 

n c (v) = «^fcK-4 i )^c(t- [fe ^ +i - 11 ) 

k 

1-2 

+ N c (rf l ^ m ' l - m -V) 

+ ( K -l)(l-l)N c (^ l -V). (88) 

The recurrence relation (|88|l reflects a single contrac- 
tion step according to the rules iJSSJi-d) and gives a sum 
of terms each containing one matrix M less than the orig- 
inal trace. Repeated such contraction steps give a sum of 
an ever increasing number of terms. After L(v) — 1 steps 
every summand will be reduced to N c (if) with v[ = I, 
v[ = for / > 2, which due to (|S3|l and 1 (85) 1 just equals 
unity. Consequently, N c (v) gives the number of terms 
in the sum at the final stage, and is thus appropriately 
called "the number of contractions" . 

Remarkably, the numbers of contractions N c (v) and 
the numbers of structures N(v) are related by 

N c (v) = N(v) 1 /' . 89 

L(v) 

With this identification, the recursion relations for both 
quantities, (H, (gDJ for N(v), and JHJJ for N c (v), 

coincide. When comparing, note that we may substitute 
Vk — Ski = v ^' l ^ k+l -)- l. A constant proportionality 
factor in (|89l) was chosen to satisfy the initial condition 
N(v') = N c (rf) = 1. In view of JgJJ, the series for K(t) 
obtained from periodic-orbit theory l|29|) and the cr-model 
(|84|l agree term by term. 

VI. CONCLUSIONS AND OUTLOOK 

Within the semiclassical frame of periodic-orbit theory, 
we have studied the spectral statistics of individual fully 
chaotic (i.e. hyperbolic and ergodic) dynamics. Central 
to our work are pairs of orbits which differ only inside 
close self-encounters. These orbit pairs yield series ex- 
pansions of the spectral form factor K(t), and our series 
agree with the predictions of random-matrix theory to all 
orders in r, for all three Wigner-Dyson symmetry classes. 
Note that we do not require any averaging over ensembles 
of systems. Moreover, we find a close analogy between 
semiclassical periodic-orbit expansions and perturbative 
treatments of the nonlinear sigma model. 

Important questions about universal spectral fluctua- 
tions remain open. The perhaps most urgent challenge 
is to go beyond the range of small r, and treat r > 1. 

The precise conditions for a system to be faithful to 
random-matrix theory remain to be established. We cer- 
tainly have to demand that the contributions of all or- 
bit pairs unrelated to close self-encounters mutually can- 
cel. While one may expect such cancellation for generic 
systems, there are important exceptions. For dynam- 
ics exhibiting arithmetic chaos, strong degeneracies in 
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the periodic-orbit spectrum give rise to system-specific 
contributions to the form factor; hence the systems in 
question deviate from random-matrix theory |3(lj |. On 
the other hand, for the Sinai billiard and the Hadamard- 
Gutzwiller model, system-specific families of orbit pairs 
found in |lflj and |29j |. respectively, do not prevent uni- 
versality. In order to formulate the precise conditions 
for the BGS conjecture, one has to clarify when non- 
universal contributions may occur. 

Moreover, a better justification is needed for neglect- 
ing the difference between stability amplitudes and peri- 
ods of the partner orbits. So far, such a justification is 
only available for Sieber/Richter pairs in the Hadamard- 
Gutzwiller model US El- 

The study of "correlated" orbit pairs opens a rich va- 
riety of applications in mesoscopic physics. Recent re- 
sults concern matrix-element fluctuations [Tsj and trans- 
port properties such as conductance, shot noise, or de- 
lay times |37|. In the latter cases, the relevant trajecto- 
ries are no longer periodic, and even e.g. quadruples of 
trajectories have interesting interpretations. While pre- 
vious work was restricted to the lowest orders in series 
expansions of the quantities in question, our machinery 
of encounters and permutations, together with intuition 
drawn from field theory should allow to attack the full 
expansion. 

Finally, one might wish to go beyond Wigner's and 
Dyson's "threefold way" and extend the present results to 
the new symmetry classes |38| , of experimental relevance 
for normal- metal/superconductor heterostructures; first 
steps are taken in [3!j . Further possible applications con- 
cern localization, a clarification of open problems in the 
nonlinear sigma model |4(j. and the crossover between 
universality classes ^| l4l| . 
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SFB/TR12 of the Deutsche Forschungsgemeinschaft is 
gratefully acknowledged. We have enjoyed fruitful dis- 
cussions with Gerhard Knieper, Jurgen Miiller, Dmitry 
Savin, Martin Sieber, Hans-Jiirgen Sommers, Dominique 
Spehner, and Martin Zirnbauer. 



end of the encounter, and only consider the stable and 
unstable separations s°, u° therein. For homogeneously 
hyperbolic dynamics, i.e. A(x, t) = e At for all x and 
t, the separations inside V c are given by s° = Sje _At ", 
!ij = Uje xtu with t u denoting the time difference between 
V and V c . 

We recall that the encounter ends when the first of 
the unstable components, say the Jth one, reaches ±c 
such that u c j = uje xtu = ±c. All I — 1 possibilities 
J = 1,2, ... ,1 — 1 and the two possibilities for the sign 
u e j/c = ±1 give additive contributions Ij to the integral 
I|A1(I . Each of them is easily evaluated after transforming 
the integration variables from Sj, Uj to s°, Uj (with j ^ 
J), s c j, and t u = j In -p^yy- The Jacobian of that transfor- 
mation equals Ac. The new coordinates are restricted to 
the ranges {— c < s° < c, — c < u° < c, for j ^ J}, —c < 
s c j < c, < t u < t cnc , and determine the action differ- 
ence as AS = J2 3 s j u j = Ej^j s j m j =t s j c - We thus 
obtain 




- A(27rft)'- 2 2ftsin^- ; (A2) 

note that the divisor t enc was canceled by the i„-integral; 
moreover, the 2(1 — 2) integrals over s°,u°, of the form 
already encountered in (|25|l . gave the factor (2irh) l ~ 2 . 

2 

Most importantly, the factor sin j- , provided by the in- 
tegral over Sj, is a rapidly oscillating function of c and h, 
annulled by averaging over these quantities; such rapidly 
oscillating terms are essentially spurious and would not 
appear if smooth encounter cut-offs were used (instead 
of our |s| < c,\u\ < c). At any rate, the integral (|A1|> . 
just the 2(1 ~ l)-fold of (|A"2)l . vanishes as h — > 0. 



APPENDIX A: INTEGRALS INVOLVING l/t enc 

We want to evaluate the integral 

f d^sd^u ^—c iAS / h (Al) 

for an Z-encounter. The integration goes over the 2(1 — 1) 
stable and unstable coordinates Sj, Uj. These variables 
determine both the duration t cnc (s,u) of the encounter 
in question and its contribution to the action difference 
AS* = J2j s j u j- We shall show that the integral may be 
neglected in the semiclassical limit. 

The key is the following change of picture: So far, all 
Poincare sections V inside a given encounter were inte- 
grated over; we thus had to divide out the duration t enc . 
Instead, we may single out a section V c , fixed at the 



APPENDIX B: EXTENSION TO GENERAL 
HYPERBOLICITY AND / > 2 

So far, we mostly restricted ourselves to two- 
dimensional homogeneously hyperbolic systems. To gen- 
eralize, we reason similarly to Ref. [ji|> where only 
Sieber/Richter pairs were considered. 

1. General hyperbolicity 

First, we shed the restriction to "homogeneously hy- 
perbolic" dynamics, for which all phase space points x 
have the same Lyapunov exponent A and stretching fac- 
tor A(t) = e xt . We shall now lift our reasoning to general 
hyperbolicity, where the stretching factors A(x, t) do de- 
pend on x. In such systems the Lyapunov exponents of 
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almost all points still coincide with the x independent 
"Lyapunov exponent of the system" , whereas each pe- 
riodic orbit may come with its own Lyapunov exponent 
p5|. Most importantly, the divergence of the stretches 
involved in an encounter depends on the local stretching 
factor of that encounter, rather than the Lyapunov expo- 
nent of the system. Our formula (|12fl for the encounter 
duration can only be read as an approximation, and that 
approximation is now to be avoided. We will thus al- 



with 



l0W *enc( X al' S c 



to depend not only on the stable 



and unstable separations s c 



but also on the phase- 



space location of the piercing x Q i chosen as the origin 
of the respective Poincare section. The changes arising 
will be important only for showing that the contribution 
arising from the 7^ — integrals of Appendix 1X1 vanishes: 
recall that for the contributing terms all occurrences of 
t onc mutually cancel. 

When generalizing the statistics of encounters of III El 
we now have to discriminate between piercing points 
x = {xn, X21, . . . , xyi} as well. Given that encounter 
stretches are separated by non-vanishing loops, these 
piercing points are uncorrelated. The analogue of wt 
will thus be a density w. r. t. x, s and u, differing from 
H17J1 only by t enc being a function of x, and by a normal- 
ization factor jjy. 

Preparing for a careful average over periodic orbits we 
first introduce the density p 7 (x, s, u, t) of piercing points, 
separations and piercing times of one fixed orbit 7, 



'(x, s, u,t) = JJ 5($ tal (zo) - Xai) 



a=l 



in 



x II S (*t^( z o) _ x "i - S aj e s (x al ) - u a je u (x a i)) . 

3=2 

Here, Zq denotes an arbitrary point of reference on 7 
and <f> t (zo) is the image of Zq under evolution over the 
time t; if the j-th stretch of the a-th encounter is al- 
most time-reversed w.r.t. the first one, we have to re- 
place $t(zo) — * T$t(zo). In analogy to lll El we integrate 
over the piercing times and divide out the encounter du- 
rations, obtaining a density of piercings and stable and 
unstable separations only, i.e. 



w 7 (x, s, u) 



J d L t p->(x,s,u,t) 

Oct ^onc( X ctlj s ai u a) 



The time integrals can be split into one over < in < T, 
and integrals over the differences t' a j = t a j — t\i of all 
other piercing times from the first one, the latter with 
the same minimal distances as in III El Using the group 
property $ taj (z ) = $ tll ($ t / (z )), we may thus repre- 
sent w 1 as the average of an observable f(z) along 7, 



'(x,a,u) = i J *n/($ tll (z ))= [/], 



/(*) = 



r v 

d L -H'l[5($ t , ai (z)-x al ) (Bl) 

^ a=l 



3=2 

The periodic-orbit sum for the form factor may now 
be written as (compare J22I) 

K /v^ N(v) 

\ V 

x / dV(x) / d L - v s d L - v u e iAS / h 



^|^ 7 | 2 5(T-T 7 )[/]. 



where the x-integral is over V points x Q i in the energy 
shell, i.e. cF/i(x) = LlLi d 4 x al 5(H(x al ) - E). 

We have not used the sum rule of Hannay and Ozo- 
rio de Almeida, except for the diagonal part. Instead, 
we invoke the equidistribution theorem |26j frecall lTTA^I . 
which says that ensembles of periodic orbits, weighted 
with their stability, behave ergodically. More precisely, if 
we average an observable /(z) (i) along a periodic orbit 
7 and subsequently (ii) over the ensemble of all such 7 
(inside a small time window and weighted with |^4 7 | 2 ), 
we obtain an energy-shell average, 



V|^ 7 | 2 5(T-T 7 )[/]. 



d/u,(z) 

n 



I AT 



/(z) = m . 

(B2) 

For the observable given in (jBlfl . the energy-shell av- 
erage can be evaluated provided the dynamics is mixing 
PH . i.e. if for two observables g(z), h(z) we have 



km g(z)h($ t (z)) 

t — >-oo 



gh. 



(B3) 



Physically, l|B3|l implies that for sufficiently large times 
t, we may neglect any classical correlations between z 
and its time evolved ^((z), and hence replace <E>t(z) by 
a phase-space point z' and average over all z'. We can 
then disregard correlations between subsequent piercing 
points with time differences at least of the order Ehren- 
fest time. Using (|B2|) . repeatedly invoking 1|B3(I for the 
product of S- functions in IjBljl . and subsequently inte- 
grating over t' a j as in lUFI we obtain 



^|^ 7 | 2 <5(T-T 7 )[/], 



'AT 



L-l 



(L-1)\UA^ L 



(B4) 
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which as expected coincides with wt(s, u) of Eq. I|17fl . up 
to division by f2 y and the x-dependence of t" nc (x, s , u )- 
The ionc-independent terms in the multinomial expan- 
sion of (|B4|) yield the same contribution to the form factor 
as before, since the divisor Q v is canceled by integration 
over x. All other contributions can be neglected in the 
semiclassical limit: they are either of a too low order in 
T or proportional to integrals of the form 



c£/i(x) 



i l - x sd l - x 



1 



^enc (^i It) 



JAS/H 



which we reveal as negligible similarly as ijAlfl . For each 



contribution Jj , we transform from x, s, u to phase-space 
points x° and separations s°, u° (u c j = ±c fixed) inside a 
Poincare section V c in the encounter end, and the sepa- 
ration t u between V and V e . For general hyperbolic dy- 
namics, the stable and unstable coordinates are related 
by s° = A(x,t u )~ 1 Sj and = A(x,t u )uj] see (Q). The 
Jacobian 7 now reads x( x °)c with the local stretching rate 
- rfln| ^ (x ' t)l ~M- We thus obtain 



defined as x{&t{*-)) 



n ( f >i*/i<i^ 1 '■ 



1 



^cnc( x °! S c , U c ) Jq 



t enc (x e ,s c ,-(i c ) 



dt u , (B5) 



coinciding with (|A2|) since the energy-shell average of the 
local stretching rate yields the Lyapunov exponent of the 
system A pEJ ]. 

The Jacobian x( xC )c has an interesting physical in- 
terpretation [l^. If we shift our Poincare section along 
the orbit, the piercing points travel, changing their un- 



stable coordinates with the velocity 



duj 
ill 



Fig. E| The Jacobian x(x e )c thus gives the velocity in 
the end of the encounter region. Restricting ourselves to 
V c and multiplying with the above velocity, we simply 
measure the flux of piercings through the line uj = ±c. 
Since each point has to traverse that line, our transforma- 
tion indeed provides an alternative counting of piercings. 
Note that the unstable coordinates temporarily shrink 
rather than grow if the local stretching rate is negative, 
and thus may traverse the line Uj — ±c several times. 
Due to the asymptotic growth of \uj\, there is one more 
traversal with growing \uj\ (and positive contribution to 
(|B5ll than with shrinking \uj\ (and negative contribu- 
tion); hence only one contribution remains effective. 



-A(x,i tl )- 1 ^J§2M«ll u e 



7 In particular, we have — —^ y ^,^ UJ u,j 

— A(x, t u )~ 1 x(x e )u e j with u e j = ±c, where we used that 
<E>t u (x) = x c . The factor A(x, t u ) _1 is compensated by the re- 
maining transformations uj — > * J) and Sj — » 




FIG. 13: Motion of piercing points through Poincare section 
V of a 2-encounter. As V is shifted, unstable components 
grow and stable ones shrink, traveling on a hyperbola AS = 
su; arrows denote the direction of motion. At end of the 
encounter, piercing points traverse the line u = c. Negative 
s, u are not shown. 



2. More than two freedoms 

For dynamics with any number / > 2 of degrees of 
freedom, the Poincare section V at point x is spanned by 
/ — 1 pairs of stable and unstable directions e|(x), ej^(x) 
(k = 1,2, ...,/— 1). A displacement <rx inside V may 
thus be decomposed as 



/-i 

£ 

fe=i 



(**e|( x )+«* e k( x )) 



compare JSJ). Each pair of directions comes with separate 
stretching factors Afc(x, t) and Lyapunov exponent Afc. 
The directions are mutually normalized as |18| | 



e«(x) Aef(x) 



<(x)Ae;(x) = e»(x)Ae i s (x) = 0, 

where (x) A (x) = 1 is a useful convention, whereas 
all other relations follow from hyperbolicity. 

Writing out the additional index fc, the uniform pier- 
cing probability (see III A|) reads ^ Yll~l dskdiik- The 
encounters (defined by \sjk\, \ujk\ < c, for all j, k) have 
heads and tails with durations 



, 1 , c 
t u = mm <^ — In — — , 

3,fc [Afc \Ujk\ 



, 1 , c 
t„ = mm \ — In — — - 

j,k I Afc Sofc 



(cp. I12|l and contribute to the action difference an 
amount given by AS — J2jk s jk u jk, with Sjk, Ujk defined 
similarly to lll C 2l The integral over s a jk , u a jk in the sec- 
ond line of (J2U yields (2Trh) < - L ~ v ^ f ~ 1 \ which is just what 
we need since Heisenberg time now reads Th = H^fkr=i ■ 
Given that the encounter ends as soon as one unstable 
component ujk reaches ±c, the -r^ — integral of Appen- 
dices ^ and IB II is split into components Ijjr, with A 
replaced by X K , and x(x) by Xif ( x )- 
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APPENDIX C: ACTION CORRELATIONS 

The semiclassical form factor |(SJ can be written in 
terms of an "action correlation function" Q , 



P(y) 



7,7' 



x 5 [y 



V. 



K(t) = t I P(y) e iy/T dy , r > 



(CI) 



Using the density of action differences P^(AS) l|19|) . we 
have evaluated the contributions to P(y) which arise from 
diagonal pairs and orbits (7,7') differing by reconnec- 
tions in close self-encounters. Collecting the terms rele- 
vant for the form factor we obtain 



P(y) 




iny + 1 £in2y m | I 

\v\ ft y ,tn 



unitary 
orthog. 



(C2) 



Random-matrix theory predicts, through the inverse 
Fourier transform of the above (|C1(I . 



Pgue(v) 
Pgoe{v) 



v 



2%) 



4 /sin(y/2) 



A ( o 

H 1 cos y si y — cos 2y si 2y 

Try \ 

+ cos y sin y [2 Ci(2y) - Ci y] ) 



with siy and Ci y the integral sine and cosine |42fl . 
There appears to be, on first glance, a contradiction 
between periodic-orbit theory and RMT. However, for 
both symmetry classes the respective results differ by 
smooth functions of the real variable y, smooth imply- 
ing continuous derivatives of all orders. According to the 
Riemann-Lebesgue theorem, the respective Fourier trans- 
forms (jCl I) have identical r-expansions of K(t). The 
even and odd parts of that expansion are respectively 
caused by the logarithmic and modulus terms in l|C2fl . 



APPENDIX D: ENCOUNTER OVERLAP 

So far, we have confined ourselves to encounters whose 
stretches are separated by intervening loops, i.e. do not 
overlap. To justify this, we shall show that encounters 
without intervening loops do not contribute to the form 
factor. The overlap of stretches involved in different en- 
counters has already been treated in 111 El We will show 
that overlap of antiparallel stretches (if not prohibited 
for dynamical reasons, as in the Hadamard-Gutzwiller 




FIG. 14: Encounters with almost self-retracing reflections: a) 
orbit scheme in configuration space; b) example for cardioid 
billiard. 



model) can be regarded as the reflection of a single stretch 
at a hard wall, e.g. in billiards. 

The case of overlapping parallel stretches is the most 
complicated one. As shown in III El two subsequent par- 
allel stretches are separated by a non-vanishing loop if 
the time £12 between the respective piercings fulfills 



£12 > £c 



t s -\- t u . 



(Dl) 



On the other hand, we shall see that any encounter with 
ti2 < t enc is accompanied by several other encounters, all 
leading to the same partner orbit; two of these encoun- 
ters obey (|D1(I . This means that blindly considering all 
encounters and associating with each of them a partner 
orbit we would count certain orbit pairs several times! It 
also follows that imposing the condition 1|D1|) we do not 
loose any orbit pair. Moreover, to count each orbit pair 
once, we need a more restrictive condition than (|D1|) . ft 
will be shown that this condition is of the form 



tl2 > tene + ifr 



(D2) 



where ^fringe is a certain "fringe duration" . It turns out 
that this boundary leads to the same value of the form 
factor as IjDlfl . which is why the simpler condition I|D1|) 
was used in the main part. 



1. Overlap of antiparallel encounter stretches 

First, let us consider two almost mutually time- 
reversed encounter stretches not separated by an inter- 
vening loop. Such a scenario is only possible if the orbit 
undergoes a nearly self-retracing reflection from a hard 
wall. After the reflection, the particle will for some time 
travel close to the pre-collision trajectory, such that tech- 
nically an antiparallel 2-encounter is formed (Fig. I14f> . 
however with just one loop and two "ports" attached. 
As shown in [l4| , no partner can be connected to such an 
encounter; formally attempting to construct a "partner" 
one obtains either the original orbit or its time-reversed. 

The proof becomes surprisingly simple if we use sym- 
bolic dynamics. Here, periodic orbits are fixed uniquely 
by sequences of symbols, e.g. denoting in certain bil- 
liards the pieces of the boundary where the orbit is re- 
flected (see e.g. ;43j); symbol sequences of periodic orbits 
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are defined modulo cyclic permutations. Even loops or 
encounter stretches can be assigned a sequence of sym- 
bols, which remains unchanged if the loop or stretch is 
slightly deflected. Given a Sieber/Richter pair, the orbit 
7 must have a sequence 7 = £\££2& , £\ and £2 denot- 
ing the two loops and £ and £ two almost time-reversed 
encounter stretches 0, |44| . The symbol sequence £ is 
obtained from £ by "time reversal" ; in the above example 
we simply have to revert the order of symbols. The part- 
ner has one loop inverted in time, and is thus described 
by the symbol sequence 7' = C\£Li£. 

As announced, both the symbol sequences of the two 
partner orbits turn out equal if the loop with symbol se- 
quence £2 is absent. A generalization to systems without 
symbolic dynamics is given in [Til ITHl flit Il8l | . 

We thus see that almost self-retracing orbit pieces as 
in Fig. 1141 each have to be regarded as one single stretch, 
folded back into itself. Clearly this carries over if they 
form part of a larger Z-encounter. 



2. Overlap of parallel encounter stretches 

Now turning to the parallel case, we first would like 
to establish an intuitive picture of "overlapping" parallel 
encounter stretches. 

For that purpose it is worthwhile to momentarily come 
back to parallel 2-encounters (even though these give rise 
to pseudo-orbits as partner orbits). The close parallel 2- 
encounters considered previously have two external loops 
attached. We shall here get to know a qualitatively dif- 
ferent kind of "parallel 2-encounter" with but a single ex- 
ternal loop. Correspondingly, there is only a single intra- 
encounter orbit stretch which remains "close to itself" 
until it exits into the external loop. Such a 2-encounter 
looks like an n-fold revolution close to some periodic or- 
bit 7; the simplest case is n — 2 but multiple revolutions 
with any integer n are possible (see Fig. II 51 for n = 3). 

Long periodic orbits with such close self-encounters are 
straightforward to construct, at least for systems with 
symbolic dynamics. To see that, we employ Fig. 116b 
which represents n = 2. Two Poincare sections mark 
beginning and end of the encounter where the orbit 7 
in question respectively enters from and departs into the 
external loop. On the other hand, these sections cut the 
encounter into two parts A and B; a bundle of three mu- 
tually close substretches traverses A, while two appear in 
B. The sequence of traversals of the two parts is ABABA. 
If the underlying system has symbolic dynamics, A and 
B each stand for a string of symbols and the foregoing 
sequence ABABA becomes a symbolic representation of 
the full encounter stretch while the small periodic orbit 
7 by itself has the code 7 = AB. With £ the symbol 
string for the external loop, the whole of 7 is represented 
by 7 = ABABAC. We may imagine that orbit to have 
arisen from another one with two external loops (sym- 
bol sequence ££'££ with £ for two parallel encounter 
stretches and £',£ for the external loops), by shrink- 




FIG. 15: Parallel 2-encounter with internal triple revolution 
about periodic orbit 7 



ing away the loop £' and thus letting the two encounter 
stretches merge into one. The special structure £ = ABA 
lets the merging stretches overlap in A and entails the or- 
bit 7 = ABABA£ in question, and that special structure 
obviously restricts the length of the overlap A to less than 
half the length of £. Otherwise, if the "overlap" extends 
over more than half the length of £ we have two pieces 
£ = (AB) n ~ 1 A overlapping in (AB) n ~ 2 A and merging 
to (AB) n A, now n > 2, which latter then represents a 
close 2-encounter with n revolutions about the shorter 
periodic orbit 7 = AB. 

Sticking to the double-revolution case of Fig. Ilfih , one 
might wonder why it is at all appropriate to speak of 
a 2-encounter, given the closeness of three substretches 
in part A. One justification was in fact already men- 
tioned: We may imagine the encounter in question gen- 
erated by shrinking away one external loop from a stan- 
dard 2-encounter. Another motivation is that of the three 
close substretches in part A (n + 1 in the case of n revo- 
lutions about 7) only two are independent; the third (the 
n — 1 others) are uniquely determined by two, calculable 
through the linearized equations of motion. 

Clearly, if one tries to reconnect (sub-) stretches of a 
parallel 2-encounter one does not get an acceptable part- 
ner orbit of 7 but a pseudo-orbit. However, the foregoing 
discussion carries over to parallel stretches in arbitrary 
Z-encounters. For instance, by shrinking away an exter- 
nal loop from a parallel 3-encounter 3 we arrive at Fig. 
116b . Here the full line revolves twice about some 7 and 
thus is the "merger" of two "overlapping" parallel en- 
counter stretches as in Fig. I16fa . while the dashed line 
revolves only once around 7. The 3-encounter of 7 in 
question thus includes a "merged" stretch ABABA and 
an "isolated" stretch £ = ABA. 

Now a partner orbit 7' may arise wherein the dashed 
stretch has pinched one revolution about 7 from the full 
one (see Fig. \Wt ) , a fact we now turn to explain. Recall 
that we can construct a partner orbit by changing connec- 
tions inside the encounter. To do so, we place a Poincare 
section within the encounter, say in region B as in Fig. 
116b . The orbit 7 will pierce through that section in three 
points; the whole of 7, the external loops included, is thus 
divided into three parts. To determine the partner 7', 
we change the connections between these parts, and thus 
their ordering. The change of connections can be inter- 
preted as cutting out one of the parts and reinserting it 



2G 



a) 




B 

b) . \ c) 




B B 



FIG. 16: Overlap ol parallel encounter stretches: a) An orbit 
7 comes close to a shorter orbit 7 (dotted line); the "merged 
stretch" (full line) approaches 7 three times in the region A 
and twice in the region B. b) An "isolated stretch" (dashed 
line) of 7 approaches A twice and B once. Also depicted is a 
Poincare section V inside the region B. c) The partner 7' has 
one repetition of 7 shifted between the merged stretch and the 
isolated stretch, d), e) Two shorter 3-encounters (depicted by 
thick lines) can be found in the region A. They pierce three 
times through a Poincare section V inside A. 



between the two other parts. In the case depicted in Fig. 
ITfib . the merged stretch pierces through V twice, and the 
part of 7 between these piercings is one revolution about 
7. The third piercing belongs to the isolated stretch. The 
partner 7' is thus obtained by cutting out one revolution 
about 7 from the merged stretches and inserting it inside 
the isolated one; see Fig. iPol . To translate to symbolic 
dynamics, we label the two remaining loops of 7 by C\ 
and £ 2 such that 7 = d (ABABA)C 2 (ABA). The part- 
ner orbit 7' has one 7 = AB shifted between the two 
stretches, i.e. 7' = C 1 (ABA)C 2 (ABABA). 

It is important to note that the encounter shown in 
Fig. 116b may be represented by three pieces A, two from 
the merged stretch and one from the isolated stretch. 
In other words, it suffices to consider a shorter encounter 
made of three independent traversals of the region A. For 
example, we may take an encounter consisting of the first 
and second traversal of A by the merged stretch and the 
latest traversal by the isolated stretch; these traversals 
are depicted in Fig. I16H by thick full and dotted lines. 



We shall refer to the full 3-encounter and its shorter rep- 
resentative as to the "long" and "short" one. 

The choice of the three representative traversals in the 
short encounter is not arbitrary First, we have to make 
sure that the encounter is limited by the boundaries of 
A. For the above example, we easily check that outside 
these boundaries, at least one of the traversing stretches 
starts to depart from the others. We have thus identified 
a shorter 3-encounter whose stretches do not overlap. 

Second, we want our short encounter to give rise to 
the same partner orbit as the long one. To check this, 
let us place a Poincare section somewhere inside A; see 
Fig. I16H . Again, the two piercings of the merged stretch 
are separated by one revolution of 7, and there is one 
piercing by the isolated stretch. Just like for the larger 
encounter, the partner has one revolution of 7 cut out of 
the merged stretch and reinserted into the isolated one. 

Both restrictions are fulfilled only for one further short 
3-encounter, consisting of the second and third traversal 
by the merged stretch and the first traversal by the iso- 
lated stretch; see Fig. 116b . 

For a first generalization we drop the above restriction 
on the length of the "overlapping part" within the merged 
stretch of a parallel 3-encounter. In general, all three 
encounter stretches must have the form £ = (AB) n ~ 1 A, 
n > 2, where the overlapping part is given by (AB) n ~ 2 A 
(the above case referred to n = 2). The merger of two 
of these reads (AB) n A and executes n revolutions about 
7 = AB, whereas the isolated stretch remains as £ = 
(AB) n ~ 1 A and revolves only n— 1 times. Next, we would 
like to generalize to parallel Z-encounters with larger If 
two of the stretches involved merge into one, this leaves 
I — 2 isolated stretches with symbol sequences as above. 
As before, the partner orbit 7' has one revolution about 
7 = AB cut out of the merged stretch and inserted into 
one of the isolated stretches. In addition, the ordering 
of the intervening loops may be changed. Finally, given 
time-reversal invariance, the isolated stretches may also 
have symbol sequences time-reversed w.r.t. (AB) n ~ 1 A. 

For all three generalizations just mentioned, we can 
identify a shorter encounter without overlap, consisting 
of I traversals of the region A, and that shorter represen- 
tative suffices to get all partner orbits 7'. To do so, we 
divide the stretches of the larger encounter in two groups. 
For each stretch of the first group, the earliest traversal 
of A will be included in the shorter encounter. For each 
stretch of the second group, we choose the latest traversal 
of A. For the "merged" stretch also the second traver- 
sal of A needs to be included (or the next-to-last, if the 
merged stretch forms part of the second group) . These I 
traversals were chosen such that before the beginning of 
A, the traversals of the first group deviate from the oth- 
ers, while the traversals of the second group do so after 
the end of A. One further short encounter is built out of 
the latest traversals by the first group, and the earliest 
traversals by the second group. As before, we disregard 
the long encounter in favor of the shorter ones. Note that 
the division of stretches in two groups is not arbitrary; 
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a) 




FIG. 17: 3-encounter (in the box) with a "fringe" where only 
two antiparallel stretches remain close. Stretches are sepa- 
rated by a loop in a) but not in b); an almost self-retracing 
reflection arises in the latter case. A partner orbit (dashed 
line) is obtained only in case a). 



the groups must be chosen to guarantee that the shorter 
encounters give rise to the same partner as the larger 
one. (See jH3 for details, and jH Hi for an alternative 
description of overlap in parallel 3-encounters.) 

To sum up, each orbit pair generated by overlapping 
parallel stretches within some encounter can be thought 
of as originating from any of two shorter encounters with- 
out overlapping stretches. One of these two shorter en- 
counters will, in the following subsection, be excluded 
and a one-to-one correspondence between orbit pairs and 
encounters will be reestablished, in effect by the condi- 
tion 1D2I) for the separation of piercings. As regards the 
form factor that condition will turn out semiclassically 
equivalent to the weaker condition I|D1(I . i.e. the require- 
ment of non-vanishing loops we had worked with in the 
main part of this paper. 



3. Fringes 

Recall that we define an Z-encounter as a region inside 
a periodic orbit in which I orbit stretches come close up 
to time reversal. Attached to the sides of the encounter 
are "fringes" in which only some of the I stretches remain 
close while others have already gone astray, as shown in 
Fig. 117b . We shall now show that these fringes do not 
affect the spectral form factor. 

As a first example, let us assume that two antiparallel 
stretches remain close after the end of the encounter, as 
in Fig. I17h . If there is no intervening loop (see Fig. 117b). 
the orbit has to undergo an almost self-retracing reflec- 
tion, just like the one studied in Appendix ID II Since no 
connections can be switched between the two stretches 
involved we have to disregard such encounters. 

On the other hand, if the "fringe" stretches are sepa- 
rated by an external loop, a partner orbit is obtained as 
usual, by reshuffling connections inside the encounter, see 
Fig. 117b . To make sure that there is an intervening loop, 
we impose minimal separations between piercing points, 
similar to those used for excluding overlap between en- 
counters in lll El After two antiparallel stretches (labelled 
by j and j + 1) have pierced through a Poincare section 
V with unstable coordinates u. 
remain close for a time 4- In 



and u Q j+i, they will 
This duration 



contains both the time t u l|12|) till the end of the a-th en- 
counter, and an additional time span after the end of the 
encounter, which gives the duration of the fringe. Using 
the unstable coordinates of piercing through a sec- 
tion in the end of the encounter (depending on u a j via 
u c a j = u a je u ), this additional time may be written as 

ihlT^ % r. 

Hence, the minimal time difference 2t u between pierc- 
ings demanded in III El has to be incremented by a time 
fringe = T In I -e _% r, depending on u (or, equiva- 

lently, the coordinates u defined in III C 2|l only via u e 
(or, equivalently, u e ). Similarly, the minimal separations 
related to stable coordinates have to be incremented by 
an amount purely depending on the stable coordinates s° 
in the beginning of the encounter. The contribution of 
each encounter to the total sum of minimal separations 
^cxci can therefore be written in the form 



,a 

L cxc\ 



U e " nc + At«( s b ) + At> e 



(D3) 



with At s and At u respective functions of s b and u c only 
(whose explicit form will not be needed in the following) . 
The numerator in the density of phase-space separations 
wt(s,u) H18|) has to be modified accordingly. 

By reasoning as in Subsection III Fl we see that only 
those terms of the multinomial expansion of the numera- 
tor in wt{s,u) contribute which involve a product of all 
i" xcl . They now have to be written as (compare l|23(l ) 



wf atI (s,u) 
L 



= m 



= Kv) 



L-V 



n L cx, 
I t° 



L-V 



n 



1 + 



Ae + At" 



'a'cn 



and contribute to the form factor as (compare 12411 ) 



KTh(v) 



T 

af a h\ 



L-V 



n 



■jL-V „ jL-V„ 



At°(s b ) + AtZ(u°) \ i E . Sa .„. 
1 H ; e R 3 3 

1 fa ' 



(D4) 



The integral in (|D4|) coincides with the one in (|24Jl up to 
the fringe corrections oc -rsp- and oc ^r 4 - , respectively. By 
reasoning similar to Appendix one easily shows that 
the resulting integrals vanish in the semiclassical limit. 
For the integrand js JL 1 the rapidly oscillating integral 

over Sj inside \A2\ is not modified; for -75*- we have 
to consider a surface of section in the beginning rather 
than in the end of the encounter. For antiparallel orbit 
stretches, fringe corrections to the form factor are hence 
revealed as negligible semiclassically. 

Let us now consider fringes of parallel encounter 
stretches. In the absence of intervening loops, there may 
be different encounters yielding the same partner orbit, 
like the two short encounters of Appendix ID 21 depicted 
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in Fig. I16H and e. One easily sees that these short en- 
counters do not have their fringes separated by loops. For 
instance, in Fig. I16H . the two full thick stretches remain 
close after the end of the encounter region A, thus form- 
ing a fringe. They are still close when the orbit reenters 
A, such that there is no external loop. 

In order to count the resulting partner 7' exactly once, 
we choose to take into account only the encounter in Fig. 
ITBfe . and disregard the one in Fig. ITBh . To do so, we again 
have to increase the minimal distances between piercing 
points; see (ID2|) . As above, the increment if r i nge may de- 
pend only on s b and u e , given that these quantities de- 
termine how long each two stretches will stay close inside 
the fringe region. More precisely, one can even show that 
tfringc is given by a sum of two terms depending on s b and 
U e , respectively. (For instance, for the first two piercings 
of a 3-encounter, tf rw 



trmge - I ln |u° I -y A m |S* 2 -S> 3 |' 

see p30| for the derivation and a generalization to arbi- 
trary /-encounters.) Hence, the sum of all minimal dis- 
tances will still be of the form (|D3|) : the increment will 
thus not contribute to K(t) in the semiclassical limit. 

The present treatment immediately generalizes to / > 
2, by writing out the additional index k, and to non- 
homogeneously hyperbolic systems by keeping the depen- 
dence on x. In particular, At s will depend both on s b 
and the phase-space points x b in the beginning of the 
encounter, and At u on u c and x e in the end of the en- 
counter. 
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APPENDIX E: THE NONLINEAR a-MODEL 

1. Replica trick and perturbation expansion 

We here sketch the derivation of the integral represen- 
tation (|79l) for the Wigner/Dyson ensembles. For sim- 
plicity, we concentrate on the GUE. For the GOE see 
Subsection IE 21 below. We start from the "replica repre- 
sentation" of the Green function G(z) = (z — -ff) -1 , 



trG(z) = d z tr ln(G(z) _1 ) = 
= — lim —d z det G(z 

r^O T 



d z ln(det G{z)- 



(El) 



We take the Hamiltonian H as an iV-dimensional Hermi- 
tian random matrix whose Gaussian statistics is defined 
by the two moments 



Hfj, v — 0, 
^2 



(E2) 



these imply Wigner's semicircle for the mean level density 
as p(E) = {N/k\ 2 )^\ 2 -{E/2) 2 Q{\ 2 - (E/2) 2 ) and 
thus the mean spacing ttX/N at E = 0. 

The rth power of the determinant in IjElj) can be writ- 
ten as the Gaussian integral (summation convention) 

det(G(z)) r = J dip exp(±f0 at (z - H)ip a ) ; (E3) 



here ip a = {V^}, a = 1, ■ • ■ , r are A^-dimensional com- 
plex vectors of integration variables; the measure dip cx 
11^=1 111=1 dRetp^dlrmp 1 ^ comprises a normalization con- 
stant of the form (const) 1 " ~ 1; here and in the remainder 
of this Appendix we let the symbol ~ mean equality in 
the limit r — > 0; the positive/negative sign in the expo- 
nent applies to the case Imz > O/Imz < 0. 

We proceed to the retarded and advanced Green func- 
tions tr G±(E) = tr G(E ± iS) with real E and 6 J. 0. 
Their ensemble averaged product 



C(e) = (trG + (E + e/2))(trG-{E-e/2)) (E4) 

leads to the spectral correlator and the form factor ac- 
cording to H77|) . Applying the replica trick i|El(l and 
the integral representation (|E3|I twice, i.e. both for G+ 
and G_ , we get to a bosonic replica variant of the sigma 
model: the correlation function G(e) appears as a two- 
fold derivative of a generating function, 



lim 

r— >o r 



d 2 Z(e) 



C(e) 

Z{e) = / dip exp(i^tA(e + E — H)ip) 



c 
A 



de+de. 



diag (e+ + iS, e_ - iS) , 6 J, , 
diag(l,-l). 



(E5) 
(E6) 



Compared to i|E3() the number of integrations is doubled 
in the generating function Z: a doublet ip = (ip 1 ,^ 2 ) 
comprises riV-dimensional vectors ip 1 and ip 2 needed to 
represent the retarded and the advanced Green functions 
by Gaussian integrals. The diagonal matrices e, A act 
in the newly introduced advanced/retarded space like 



A 



= (-^2), arid like unity in the replica space. 
Note that A and the imaginary part of e together se- 
cure convergence of the Gaussian integrals in Z as did 
the plus/minus sign in l|E3|l . The replica index a is sup- 
pressed in the compact notation for Z. 

The GUE average is now easily done in 1E6J) since the 
random Hamiltonian appears linearly in the exponent, 

Z(e) =Jd?P expfo t A(e + E)i> - A^tr A 2 ) 

with the 2r x 2r matrix A = {A afs = V^V^A^}; 
the index a = (p, a) comprises the replica index a — 
1 . . . r and an index p — 1, 2 discriminating between "ad- 
vanced" and "retarded". We next decouple the term 
oc tr A 2 quartic in the integration variables ip by a 
Hubbard-Stratonovich transformation. To this end, we 
multiply Z by the unit-normalized Gaussian integral 
1 = AffdQ exp (f trQ 2 ), where Q = {Q af3 } is a 2r- 
dimensional anti-Hermitian matrix and J dQ denotes in- 
tegration over all its independent matrix elements; the 
latter comprise (2r) 2 independent real and imaginary 
parts such that the normalization constant is of the form 
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N = (const) r ~ 1. The shift Q 



XA/N leads to that tr(Q 2 ) = tr(A 2 



and that tr(Q s e) 



Z(e) ~Jdip JdQexp(§ tr Q 2 + A{e + E + iXQ)^) 
~JdQ exp (f tr Q 2 - N tr ln(e + E + iAQ)) 

where in the second step we have performed the Gaussian 
integration over ip. The large parameter N in the fore- 
going exponential allows for a stationary-phase approxi- 
mation of the Q-integral. Variation of the exponent, the 
so-called action, w.r.t. the independent matrix elements 
of Q yields the saddle-point equation 



Q 



iA 



e + E + iAQ 



(E7) 



In order to shorten the saddle-point analysis we confine 
ourselves to £ = 0; the final result is valid through- 
out Wigner's semicircle. The saddle-point equation then 
invites solution to zeroth order in e since the scale for 
the energy offset variables is the mean level spacing, 
e± = 0(N^ 1 ); we can even dispatch the infinitesimal 
imaginary part Im(e) = iA<5 which is only needed to 
ensure absolute convergence of all integrals. A sim- 
ple matrix-diagonal saddle point thus arises, Q ~ A. 
(Strictly speaking, there are 2 2r diagonal solutions with 
diagonal entries ±1. However, the solutions different 
from A = diag(l,— 1), have to be discarded since they 
cannot be reached without crossing the cuts of the loga- 
rithm in the action |45j). 

A continuous manifold of more general (non-diagonal 
yet compatible with the cut structure of the logarithm) 
solutions of the saddle-point equation can now be ob- 
tained by the conjugation A — * TAT -1 = Q s , where 
T G U(r, r) is a (2r)-dimensional pseudo-unitary matrix. 
(The pseudounitarity condition, TKT^ = A, is required 
to make the subse quen t integration over all configura- 
tions Q convergent |4G|.) This observation identifies the 
non-compact symmetric space U(r, r)/U(r) x U(r) as the 
saddle-point manifold. (Transformations T £ U(r)xU(r) 
commute with A and, therefore, do not affect the diag- 
onal saddle.) It may be a noteworthy property of the 
saddle-point manifold in question that the matrices Q 
thereon are no longer anti-Hermitian |46| . 

We next substitute the saddle configurations Q s = 
TAT -1 into the action and expand to order e, 

Z(e) ~ J dQ s exp(-iVtr [e(E + iXQs)- 1 ]) 
dQ s exp I — tr [eQ s \ 
dQ s exp '— tr [AQ S 



dQ s exp ( — trAQ E 



(E8) 



where now J dQ s demands integration over the saddle- 
point manifold 45] ; in the foregoing calculation we used 



(e + + e_)tr(Q s )/2 + (e+ - e_ + 2i<5)tr(Q s A)/2 = (e+ - 
e_ + 2i<5)tr(Q s A)/2; finally, in the last step we introduced 
the scaled offset s + — (e + — e_ + 2i<5)7rp. In order to 
evaluate the matrix integral above we need an explicit 
parametrization of the Q-matrices on the saddle-point 
manifold. Taylor made to perturbative calculations is 
the so-called rational parametrization 8 



? s = (1 - W)A(1 - W)- 1 ; W = 



1 {OB 



7+ V Bt o 



(E9) 



where the scaling factor l/vs + makes for notational 
convenience. The main merit of that representation is 
that the invariant measure dQ s is effectively (i.e. in 
view of the eventual limit r — > 0) the flat measure 
(s+)- r2 Hlj r dReBijdimBij = {s+)- r2 dB for r x r ma- 
trices B. Inserting the parametrization l|E9(l in the gener- 
ating function l|E8|l we expand as (1 — W)~ x = X^^Lo 
to get 



Z(e) ~ ( s +)- 2 e is4> / dBe {2i / K ^^ s+ ^ lt ^ BB " )l , 



(E10) 

where the factor k — 1 has been sneaked in for later use. 

We had shown in IV Bl that the foregoing i?-integral 
yields {const r2 +r 2 /2(s + )+r 3 /3(s + )+0(r 4 )} w i t h f l (s+) 
polynomials in l/s + . Therefore, nonzero contributions 
to the correlator C(e) according to 1E5|) arise only when 
both derivatives d 2 /de+de- act together on each of the 
three factors on the r. h. s. of ijElOfl . It is easy to see that 
the contribution from d 2 e ls+r /de + de- yields the discon- 
nected part p 2 such that upon simply setting e ls+r — > 1 
and invoking (|E5(1 we get (the non-oscillatory contribu- 
tions to) the connected part C conn (e) instead of C(e). 
The connected correlation function defined in (|77|l is thus 
obtained as R — ^3=iReC c „„„ = — ^Re-J^Hr(Z\ 



h Re d 3 

2 as + 



which leads us to l(79")l for the unitary class, with k = 1. 

2. Generalization to the orthogonal class 

For the GOE, Eq. (|E2(1 generalizes to 
(H^) = 0, 

A 2 

As a consequence, the averaged action will contain the 
sum of two contributions quartic in the integration 
variables, ^A^ (^JA^ + A^) = 2^ At A* l/ 'l'^A* y , 



8 Other parametrizations allow to conveniently do the Q s -integral 
exactly, to get the non-oscillatory part of R(s) in closed form 
instead of the asymptotic l/s+-series we are after. 
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where we have defined 



V2 



* = 




(Ell) 



Notice that the 4r-component vectors \& and VP are con- 
nected by the symmetry relation 



* = V 1 (Tl, 



(E12) 



where the Pauli matrix 01 acts in the two-component 
space of Eq. (|E11|) . To decouple the quartic \l/-term 
in the action, we introduce a Q-matrix subject to the 
constraint Q = (TiQ 1 o\. Noting that ip^(e + E)tp — 
^(e + E)^ and proceeding as in the unitary case we get 

Z{e) = Jd^JdQexp(^tiiQ 2 )+MA(i + E + i\Q)-<Z 

= JdQ exp (^tr(Q 2 ) - f tr ln(e + E + iAQ)) , 

where we made use of the symmetry of Q. A stationary- 
phase approach yields the diagonal saddle Q = A and 
its non-diagonal generalization Q s = TAT -1 , where 
T e 0(2r, 2r). (Compatibility with the symmetries 
of Q requires T T = a\T~ 1 (Ti, hence the restric- 
tion to the (pseudo)orthogonal group.) This identifies 
0(2r, 2r)/0(2r) x 0(2r) as the saddle-point manifold. 

As in the unitary case, we represent the Q s -matrices 
in the rational parametrization l|E9|l . Presently, time- 
reversal invariance implies T T = o\T~ l o\ and, as the 
symmetry of the integration variables B, 



= -a 1 B T a 1 



(E13) 



Simplifying the action to leading order in e and expand- 
ing in W, we recover Eq. I|79|) . now with k = 2. 



3. Contraction rules 

To compute the matrix integral perturbatively in 
an expansion in l/s + , one expands the exponenti- 
ated action in powers of traces ti[(BB^) 1 ], I > 2. 
This leads to a series of terms of the structure 
((tT[(BB^) 2 ]) V2 (tr[(BB^) 3 }) V3 ...), where v t are integers 
and the averaging over the quadratic action, (...), is de- 
fined by H81|) . Each term contributing to the series is then 



computed by Wick's theorem, i.e. as a sum over all non- 
vanishing pair contractions (|85b -d) of i3-matrices. We 
now briefly outline the derivation of the basic contrac- 
tion rules. In an index notation, the quadratic action 
reads S^[B,B^} = ii/n)Y, aa ' \B a a'\ 2 - This implies the 
prototypical contraction rule 



(B ryn i B^ f<i 



2i 



the ^-independence of which follows from the fact that 
in the orthogonal case, k = 2, only one half of the matrix 
elements B aa i are independent integration variables - 
a consequence of the time reversal relation (|E13I) . When 
expressed as a sum over independent integration variables 
only, the coefficient of the action effectively doubles, i.e. 
in either case k = 1, 2 the integration over matrix compo- 
nents obtains a factor i/2. Using this relation we obtain 



I 1 



r 



(tiBXtT&Y) = (B a0 X Pa B\ s Y Sl ) 



(Xp a Y c 



a(3) 



;tr(AF)), 



which is contraction rule (|85b ,~l. Rule (JHSJd) is proven in 
the same manner. Rule (JHSJi) for the GOE results from 



(tvBXBY) 



(tr{{a 1 X)B(Ya 1 )(a 1 Ba 1 ))) 
{{viX) a pBp 1 (Y<7 1 ) lS (<T 1 B<Ji)s a ) 



= +j i ((a 1 XU(a 1 Y T )p a ) 
= +i<tr(A(a 1 y T a 1 ))) 
= -I(tr(XYt)), 

where in the last step we used that Y contains an odd 
number of matrices B and B^ (i.e. that the time reversal 
relation (|E13|I implies a{Y T u\ — — YU) The proof of 
rule (JHHt) proceeds along the same lines. 
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